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ABSTRACT 

The European Spaee Agency’s Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 
and scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013. In February 2015, ESA and the 
Planck Collaboration released the seeond set of eosmology produets based on data from the entire Planck mission, including both temperature and 
polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main 
characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. 
The scienee products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds 
in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich 
clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and 
assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB 
lensing likelihood. Scientific results inelude cosmological parameters deriving from CMB power speetra, gravitational lensing, and cluster eounts, 
as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity. 

Key words. Cosmology: observations - Cosmic background radiation - Surveys - Space vehicles: instruments - Instrumentation: detectors 
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1, Introduction r 

The Planck satellite^ (Tauber et al. 2010; Planck Collaboration I 
2011) was launched on 14 May 2009 and observed the sky 
stably and continuously from 12 August 2009 to 23 October 
2013. Planck's scientific payload contained an array of 74 de¬ 
tectors in nine frequency bands sensitive to frequencies between 
25 and 1000 GHz, which scanned the sky with angular resolu¬ 
tion between 33' and 5'. The detectors of the Low Frequency 
Instrument (LFI; Bersanelli et al. 2010; Mennella et al. 2011) 
were pseudo-correlation radiometers, covering bands centred 
at 30, 44, and 70 GHz. The detectors of the High Frequency 
Instrument (HFI; Lamarre et al. 2010; Planck HFI Core Team 
2011a) were bolometers, covering bands centred at 100, 143, 
217, 353, 545, and 857 GHz. Planck imaged the whole sky twice 
in one year, with a combination of sensitivity, angular resolution, 
and frequency coverage never before achieved. Planck, its pay- 
load, and its performance as predicted at the time of launch are 
described in 13 papers included in a special issue of Astronomy 
& Astrophysics (Volume 520). 

The main objective of Planck, defined in 1995, was to mea¬ 
sure the spatial anisotropies in the temperature of the cos¬ 
mic microwave background (CMB), with an accuracy set by 
fundamental as trophy sical limits, thereby extracting essentially 
all the cosmological information embedded in the temperature 
anisotropies of the CMB. Planck was not initially designed to 
measure to high accuracy the CMB polarization anisotropies, 
which encode not only a wealth of cosmological information, 
but also provide a unique probe of the early history of the 
Universe, during the time when the first stars and galaxies 
formed. However, during its development it was significantly 
enhanced in this respect, and its polarization measurement ca¬ 
pabilities have exceeded all original expectations. Planck was 
also designed to produce a wealth of information on the proper¬ 
ties of extragalactic sources and of clusters of galaxies (via the 
Sunyaev-Zeldovich effect), and on the dust and gas in the Milky 
Way. The scientific objectives of Planck were described in de¬ 
tail in Planck Collaboration (2005). With the results presented 
here and in a series of accompanying papers, Planck has already 
achieved all of its planned science goals. 

An overview of the scientific operations of the Planck 
mission has been presented in Planck Collaboration I (2014). 
Further operational details—extending this description to the 
end of the mission—are presented in the 2015 Explanatory 
Supplement (Planck Collaboration ES 2015). This paper 
presents an overview of the main data products and scientific 
results of Planck's third release,^ based on data acquired in 
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^Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from Prance and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 

^In January of 2011, ESA and the Planck Collaboration released 
to the public a first set of scientific data, the Early Release Compact 
Source Catalogue (ERCSC), a list of unresolved and compact sources 
extracted from the first complete all-sky survey carried out by Planck 
(Planck Collaboration VII 2011). At the same time, initial scientific re¬ 
sults related to astrophysical foregrounds were published in a special 
issue of Astronomy and Astrophysics (Vol. 520, 2011). Since then, 
34 UPDATE “Intermediate” papers have been submitted for publi¬ 
cation to A&A, containing further astrophysical investigations by the 


the period 12 August 2009 to 23 October 2013, and hereafter 
referred to as the “2015 products.” 

2. Data products in the 2015 release 

The 2015 distribution of released products, which can be freely 
accessed via the Planck Legacy Archive interface (PLA),^ is 
based on all the data acquired by Planck during routine oper¬ 
ations, starting on 12 August 2009 and ending on 23 October 
2014. The distribution contains the following. 

- Cleaned and calibrated timelines of the data for each detec¬ 
tor. 

- Maps of the sky at nine frequencies (Sect. 7) in tempera¬ 
ture, and at seven frequencies (30-353 GHz) in polarization. 
Additional products serve to quantify the characteristics of 
the maps to a level adequate for the science results being 
presented, such as noise maps, masks, and instrument char¬ 
acteristics. 

- Four high-resolution maps of the CMB sky in temperature 
and polarization, and accompanying characterization prod¬ 
ucts (Sect. 8.1)"^. 

- Four high-pass-filtered maps of the CMB sky in polarization, 
and accompanying characterization products (Sect. 8.1). The 
rationale for providing these maps is explained in the follow¬ 
ing section. 

- A low-resolution CMB temperature map (Sect. 8.1) used in 
the low-7’ likelihood code, with an associated set of fore¬ 
ground temperature maps produced in the process of sepa¬ 
rating the low-resolution CMB from foregrounds, with ac¬ 
companying characterization products. 

- Maps of thermal dust and residual cosmic infrared back¬ 
ground (CIB), carbon monoxide (CO), synchrotron, free- 
free, and spinning dust temperature emission, plus maps of 
dust temperature and opacity (Sect. 9). 

- Maps of synchrotron and dust polarized emission. 

- A map of the estimated CMB lensing potential over 70 % of 
the sky. 

- A map of the Sunyaev-Zeldovich effect Compton parameter. 

- Monte Carlo chains used in determining cosmological pa¬ 
rameters from the Planck data. 

- The Second Planck Catalogue of Compact Sources 
(PCCS2, Sect. 9.1), comprising lists of compact sources 
over the entire sky at the nine Planck frequencies. 
The PCCS2 includes polarization information, and su¬ 
persedes the previous Early Release Compact Source 


Collaboration. In March of 2013, the second release of scientific data 
took place, consisting mainly of temperature maps of the whole sky; 
these products and associated scientific results are described in a spe¬ 
cial issue of A&A (Vol. 571, 2014). 

^http://pla.esac.esa.int 

^It has become the norm in CMB studies to use the COSMO 
(http://healpix.sourceforge.net/html/intronode6.htm) 
convention for polarization angles, rather than the lAU 
(Hamaker & Bregman 1996) convention, and the Planck data products 
have followed this trend. The net effect of using the COSMO conven¬ 
tion is a sign inversion on Stokes U with respect to the lAU convention. 
All Planck fits files containing polarization data include a keyword that 
emphasizes the convention used. Nonetheless, users should be keenly 
aware of this fact. 
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Catalogue (Planck Collaboration XIV 2011) and the PCCSl 
(Planck Collaboration XXVIII 2014). 

- The Second Planck Catalogue of Sunyaev-Zeldovich 
Sources (PSZ2, Sect. 9.2), comprising a list of sources de¬ 
tected by their SZ distortion of the CMB spectrum. The 
PSZ2 supersedes the previous Early Sunyaev-Zeldovich 
Catalogue (Planck Collaboration XXIX 2014) and the PSZl 
(Planck Collaboration XXIX 2014). 

- The Planck catalogue of Galactic Cold Clumps (PGCC, 
Planck Collaboration XXVIII 2015), providing a list of 
Galactic cold sources over the whole sky (see Sect. 9.3). The 
PGCC supersedes the previous Early Cold Core Catalogue 
(ECC), part of the Early Release Compact Source Catalogue 
(ERCSC, Planck Collaboration VII 201 1). 

- A full set of simulations, including Monte Carlo realizations. 

- A likelihood code and data package used for testing cosmo¬ 
logical models against the Planck data, including both the 
CMB (Sect. 8.4.1) and CMB lensing (Sect. 8.4.2). 

The first 2015 products were released in Eebruary 2015, and 
the full release will be complete by July 2015. In parallel, the 
Planck Collaboration is developing the next generation of data 
products, which will be delivered in the early part of 2016. 

2.1. The state of polarization in the Planck 2015 data 

LFI —The 2015 Planck release includes polarization data at the 
LEI frequencies 30, 44, and 70 GHz. The 70 GHz polarization 
data are used for the 2015 Planck likelihood at low-multipoles 
{i < 30). The 70 GHz map is cleaned with the 30 GHz and the 
353 GHz channels for synchrotron and dust emission, respec¬ 
tively (Planck Collaboration XIII 2015). 

Control of systematic effects is a challenging task for po¬ 
larization measurements, especially at large angular scales. We 
carry out extensive analyses of systematic effects impacting the 

2015 LEI polarization data (Planck Collaboration III 2015). Our 
approach follows two complementary paths. Eirst, we use the 
redundancy in the Planck scanning strategy to produce differ¬ 
ence maps that, in principle, contain the same sky signal (“null 
tests”). Any residuals in these maps blindly probe all non- 
common-mode systematics present in the data. Second, we use 
our knowledge of the instrument to build physical models of all 
the known relevant systematic effects. We then simulate time¬ 
lines and project them into sky maps following the LEI map¬ 
making process. We quantify the results in terms of power spec¬ 
tra and compare them to the EEP8 LEI noise model. 

Our analysis shows no evidence of systematic errors sig¬ 
nificantly affecting the 2015 LEI polarization results. On the 
other hand, our model indicates that at low multipoles the dom¬ 
inant LEI systematics (gain errors and ADC non-linearity) are 
only marginally dominated by noise and the expected signal. 
Therefore, further independent tests are being carried out and 
will be discussed in a forthcoming paper, as well as in the final 

2016 Planck release. These include polarization cross-spectra 
between the LEI 70 GHz and the HEI100 and 143 GHz maps 
(that are not part of this 2015 release - see below). Because sys¬ 
tematic effects between the two Planck instruments are expected 
to be largely uncorrelated, such cross-instrument approach may 
prove particularly effective. 

HFI —The February 2015 data release included polarization 
data at 30, 44, 70, and 353 GHz. The release of the remain¬ 
ing three HFI channels - 100, 143, and 217 GHz - was de¬ 
layed because of residual systematic errors in the polarization 


data, particularly but not exclusively at ^ < 10. The sources of 
these systematic errors were identified, but insufficiently char¬ 
acterized to support reliable scientific analyses of such things 
as the optical depth to ionization r and the isotropy and statis¬ 
tics of the polarization fluctuations. Due to an internal mixup, 
however, the unfiltered polarized sky maps ended up in PEA in¬ 
stead of the high-pass-filtered ones. This was discovered in July 
2015, and the high-pass-filtered maps at 100, 143, and 217 GHz 
were added to the PEA. The unfiltered maps have been left in 
place to avoid confusion, but warnings about their unsuitability 
for science have been added. Since February our knowledge of 
the causes of residual systematic errors and our characterization 
of the polarization maps have improved. Problems that will be 
encountered in the released 100-353 GHz maps include the fol¬ 
lowing: 

- Null tests on data splits indicate inconsistency of polariza¬ 
tion measurements on large angular scales at a level much 
larger than our instrument noise model (see Fig. 10 of 
Planck Collaboration VIII2015). The reasons for this are nu¬ 
merous and will be described in detail in a future paper. 

- While analogue-to-digital converter (ADC) nonlinearity is 
corrected to a much better level than in previous releases, 
some residual effects remain, particularly in the distortion of 
the dipole that leaks dipole power to higher signal frequen¬ 
cies. 

- Bandpass mismatches leak dust temperature to polarization, 
particularly on large angular scales. 

- While the measured beam models are improved, main beam 
mismatches cause temperature-to-polarization leakage in the 
maps (see Fig. 17 of Planck Collaboration VII2015). In pro¬ 
ducing the results given in the Planck 2015 release, we cor¬ 
rect for this at the spectrum level (Planck Collaboration XI 
2015), but the public maps contain this effect. 

The component separation work described in Sect. 9, 
Planck Collaboration IX 2015, and Planck Collaboration X 
2015 was performed on all available data, and produced un¬ 
precedented full-sky polarization maps of foreground emission 
(Figs. 21 and 23), as well as maps of polarized CMB emission. 
The polarized CMB maps, derived using four independent com¬ 
ponent separation methods, were the basis for quantitative state¬ 
ments about the level of residual polarization systematics and 
the conclusion that reliable science results could not be obtained 
from them on the largest angular scales. 

Recent improvements in mapmaking methodology that re¬ 
duce the level of residual systematic errors in the maps, espe¬ 
cially at low multipoles, will be described in a future paper. A 
more fundamental ongoing effort aimed at correcting systematic 
polarization effects in the time-ordered data will produce the fi¬ 
nal legacy Planck data, to be released in 2016. 

3. Papers accompanying the 2015 release 

The characteristics, processing, and analysis of the Planck data, 
as well as a number of scientific results, are described in a series 
of papers released with the data. The titles of the papers begin 
with ''Planck 2015 results.”, followed by the specific titles below. 

I. Overview of products and scientific results (this paper) 

II. Low Frequency Instrument data processing 

III. LFI systematic uncertainties 

IV. LFI beams and window functions 
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V. LFI calibration 

VI. LFI maps 

VII. High Frequency Instrument data processing: Time-ordered 

information and beam processing 

VIII. High Frequency Instrument data processing: Calibration 

and maps 

IX. Diffuse component separation: CMB maps 

X. Diffuse component separation: Foreground maps 

XI. CMB power spectra, likelihoods, and robustness of 

parameters 

XII. Simulations 

XIII. Cosmological parameters 

XIV. Dark energy and modified gravity 

XV. Gravitational lensing 

XVI. Isotropy and statistics of the CMB 

XVII. Constraints on primordial non-Gaussianity 

XVIII. Background geometry and topology of the Universe 

XIX. Constraints on primordial magnetic fields 

XX. Constraints on inflation 

XXL The integrated Sachs-Wolfe effect 

XXII. A map of the thermal Sunyaev-Zeldovich effect 

XXIII. The thermal Sunyaev-Zeldovich effect-cosmic infrared 
background correlation 

XXIV. Cosmology from Sunyaev-Zeldovich cluster counts 

XXV. Diffuse, low-frequency Galactic foregrounds 

XXVI. The Second Planck Catalogue of Compact Sources 

XXVII. The Second Planck Catalogue of Sunyaev-Zeldovich 
Sources 

XXVIII. The Planck Catalogue of Galactic Cold Clumps 

This paper contains an overview of the main aspects of the 
Planck project that have contributed to the 2015 release, and 
points to the papers that contain full descriptions. It proceeds 
as follows. Section 4 describes the simulations that have been 
generated to support the analysis of Planck data. Section 5 de¬ 
scribes the basic processing steps leading to the generation of 
the Planck timelines. Section 6 describes the timelines them¬ 
selves. Section 7 describes the generation of the nine Planck fre¬ 
quency maps and their characteristics. Section 8 describes the 
Planck 2015 products related to the cosmic microwave back¬ 
ground, namely the CMB maps, the lensing products, and the 
likelihood code. Section 9 describes the Planck 2015 astrophysi- 
cal products, including catalogues of compact sources and maps 
of diffuse foreground emission. Section 10 describes the main 
cosmological science results based on the 2015 CMB products. 
Section 11 describes some of the astrophysical results based on 
the 2015 data. Section 12 concludes with a summary and a look 
towards the next generation of Planck products. 

4. Simulations 

We simulate time-ordered information (TOI) for the full focal 
plane (FFP) for the nominal mission. The first five FFP real¬ 
izations were less comprehensive and were primarily used for 
validation and verification of the Planck analysis codes and 
for cross-validation of the data processing centres’ (DPCs) and 
FFP simulation pipelines. The first Planck cosmology results 
(Planck Collaboration I 2014) were supported primarily by the 
sixth FFP simulation-set, hereafter FFP6. The current results 


were supported by the next generation of simulations, FFP8, 
which is described in detail in Planck Collaboration XII (2015). 

Each FFP simulation comprises a single “fiducial” realiza¬ 
tion (CMB, astrophysical foregrounds, and noise), together with 
separate Monte Carlo (MC) realizations of the CMB and noise. 
The CMB component contains the effect of our motion with 
respect to the CMB rest frame. This induces an additive dipo¬ 
lar aberration, a frequency-dependent dipole modulation, and a 
frequency-dependent quadrupole in the CMB data. Of these ef¬ 
fects, the additive dipole and frequency-independent component 
of the quadrupole are removed (see Planck Collaboration XII 
2015 for details), the residual quadrupole, aberration, and mod¬ 
ulation effects are left in the simulations and are also left in the 
LFI and HFI data. To mimic the Planck data as closely as pos¬ 
sible, the simulations use the actual pointing, data flags, detec¬ 
tor bandpasses, beams, and noise properties of the nominal mis¬ 
sion. For the fiducial realization, maps were made of the total 
observation (CMB, foregrounds, and noise) at each frequency 
for the nominal mission period, using the Planck Sky Model 
(Delabrouille et al. 2013). In addition, maps were made of each 
component separately, of subsets of detectors at each frequency, 
and of half-ring and single Survey subsets of the data. The noise 
and CMB Monte Carlo realization-sets also included both all and 
subsets of detectors (so-called “DetSets”) at each frequency, and 
full and half-ring data sets for each detector combination. 

To check that the PR2-2015 results are not sensitive to the 
exact cosmological parameters used in FFP8, we subsequently 
generated FFP8.1, exactly matching the PR2-2015 cosmology. 

All of the FFP8 and FFP8.1 simulations are available to be 
used at NERSC (http://crd.lbl.gov/cmb-data); in addi¬ 
tion, a limited subset of the simulations are available for down¬ 
load from the PLA. 

5. Data Processing 

5.1. Timeline processing 

5.1.1. LFI 

The main changes in the LFI data processing compared to the 
earlier release (Planck Collaboration II 2014) are related to the 
way in which we take into account the beam information in the 
pipeline processing, as well as an entire overhaul of the iterative 
algorithm used to calibrate the raw data. The process starts at 
Level 1, which retrieves all the necessary information from data 
packets and auxiliary data received from the Mission Operation 
Centre, and transforms the scientiflc packets and housekeeping 
data into a form manageable by Level 2. Level 2 uses scientiflc 
and housekeeping information to: 

- build the LFI reduced instrument model (RIMO), which con¬ 
tains the main characteristics of the instrument; 

- remove ADC non-linearities and 1 Hz spikes diode by diode; 

- compute and apply the gain modulation factor to minimize 
1 If noise; 

- combine signals from the diodes with associated weights; 

- compute the appropriate detector pointing for each sample, 
based on auxiliary data and beam information, corrected by 
a model (PTCOR) built using solar distance and radiometer 
electronics box assembly (REBA) temperature information; 

- calibrate the scientiflc timelines to physical units (Kcmb), fit¬ 
ting the total CMB dipole convolved with the An beam rep¬ 
resentation, without taking into account the signature due to 
Galactic straylight; 
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- remove the solar and orbital dipole (convolved with the An 
beam) representation and the Galactic emission (convolved 
with the beam sidelobes) from the scientific calibrated time¬ 
line; 

- combine the calibrated time-ordered information (TOI) into 
aggregate products, such as maps at each frequency. 

Level 3 collects Level 2 outputs from both HFI 
(Planck Collaboration VIII 2015) and LFI and derives various 
products, such as component-separated maps of astrophysical 
foregrounds, catalogues of different classes of sources, and the 
likelihood of cosmological and astrophysical models given in 
the maps. 

5.1.2. HFI 

The HFI data processing for this release is very similar to that 
used for the 2013 release (Planck Collaboration VI 2014). The 
main improvement is carried out in the very first step in to the 
pipeline, namely the correction for ADC non-linearities. 

The HFI bolometer electronic readout, described in 
the Planck Explanatory Supplement (Planck Collaboration ES 
2015), ends with a 16-bit Analogue-to-Digital Converter. Its tol¬ 
erance on the differential non-linearity (the maximum deviation 
from one least significant bit, LSB, between two consecutive lev¬ 
els, on the whole range) is specified to be better than +1.6 LSB. 
The consequences of this feature on HFI performances had 
not been anticipated, nor did it produce any detected effect on 
ground-test data, but it proved to be a major systematic effect 
impacting the flight data. A method that reduces the ADC effect 
by more than an order of magnitude for most channels has been 
implemented. 

No changes were made to any software module involved 
in the TOI processing, from ADC-corrected TOI to clean TOI 
that are ready for qualification, calibration and mapmaking. 
However, several input parameters of the modules have been 
fine-tuned for better control of some residual systematic errors 
that were noticed in the 2013 data. 

Improvements can be assessed by comparing the noise sta- 
tionarity in the 2013 and 2015 data. Trends of the so-called total 
noise versus ring number before (black dots, 2013 release) and 
after the ADC correction (blue dots, this release) are shown in 
Fig. 1. There is a significant decrease in the relative width of 
the distribution when the ADC correction is included. For most 
bolometers, the noise stationarity is ascertained to be within the 
percent level (Planck Collaboration VIII 2015). 

For strong signals, the threshold for cosmic ray removal 
(“deglitching”) is auto-adjusted to cope with source noise, due 
to the small pointing drift during a ring. Thus, more glitches 
are left in data in the vicinity of bright sources, such as the 
Galactic centre, than are left elsewhere. To mitigate this effect 
near bright planets, the signal at the planet location is flagged 
and interpolated prior to the TOI processing. For the 2015 re¬ 
lease, this is done for Jupiter at all HFI frequency bands, for 
Saturn at v > 217 GHz and for Mars at v > 353 GHz. 

Nevertheless, for beam and calibration studies 
(see Sect. 5.2.2, Planck Collaboration VII 2015 and 
Planck Collaboration VIII 2015), the TOI of all planet cross¬ 
ings, including the planet signals, are needed at all frequencies. 
Hence, a dedicated production is done in parallel for those 
pointing periods and bolometers. In that case, in order to 
preserve the quality of the deglitching, an iterative 3-level 
deglitcher has been added in the 2015 data analysis. 


As noted in Planck Collaboration I (2014), Planck scans a 
given ring on the sky for roughly 45 minutes before moving on 
to the next ring. The data between these rings, taken while the 
spacecraft spin-axis is moving, are discarded as “unstable.” The 
data taken during the intervening “stable” periods are subjected 
to a number of statistical tests to decide whether they should 
be flagged as usable or not (Planck Collaboration VI 2014); this 
procedure continues to be used for the present data release. An 
additional selection process has been introduced to mitigate the 
effect of the 4-K lines (i.e., periodic cooler variations) on the 
data, especially the 30 Hz line signal, which is correlated across 
bolometers. It is therefore likely that the 4-K line-removal pro¬ 
cedure leaves correlated residuals on the 30 Hz line. The conse¬ 
quence of this correlation is that the angular cross-power spectra 
between different detectors can show excess power at multipoles 
around i ^ 1800. To mitigate this effect, we discard all 30 Hz 
resonant rings for the 16 bolometers between 100 and 353 GHz 
for which the median average of the 30 Hz line amplitude is 
above 10 aW. As a result, the i ^ 1800 feature has now dis¬ 
appeared. 

Figure 2 summarizes the situation, showing the fraction of 
discarded samples for each detector over the full mission. It gath¬ 
ers the flags at the sample level (blue line), which are mainly 
due to glitches (green line) plus the pointing maneuvers be¬ 
tween rings (about 8 %) and the glitch flag combination for the 
polarization-sensitive bolometers (PSBs) and secondly, at the 
ring level (black line), which are mostly due to the 4-K lines, but 
also due to Solar flares, big manoeuvres, and end-of-life calibra¬ 
tion sequences, which are common to all detectors. With respect 
to the nominal mission, presented in the 2013 papers, the main 
difference appears in Survey 5, which is somewhat disjointed, 
due to Solar flares arising from the increased Solar activity, and 
to special calibration sequences. The full cold Planck HFI mis¬ 
sion lasted 885 days, excluding the Calibration and Performance 
Verification (CPV) period of 1.5 months. Globally, for this dura¬ 
tion, the total amount of HFI data discarded amounted to 31 %, 
the majority of which came from glitch flagging. 

Details of the TOI processing are given in the 
Planck Collaboration VII (2015). 


5.2. Beams 
5.2.1. LFI beams 

As described in Planck Collaboration IV (2015), the in-flight as¬ 
sessment of the LFI main beams relied on the measurements 
performed during seven Jupiter crossings: the first four transits 
occurred in nominal scan mode (spin shift 2', l°day“^); and 
the last three scans in deep mode (spin shift 0.'5, 15'day“^). 
By stacking data from the seven scans, the main beam profiles 
are measured down to -25 dB at 30 and 44 GHz, and down to 
-30 dB at 70 GHz. Fitting the main beam shapes with an ellipti¬ 
cal Gaussian profile, we have expressed the uncertainties of the 
measured scanning beams in terms of statistical errors for the 
Gaussian parameters: ellipticity; orientation; and FWHM. In this 
release, the error on the reconstructed beam parameters is lower 
with respect to that in the 2013 release. Consequently, the error 
envelope on the window functions is lower as well. For exam¬ 
ple, the beam FWHM is determined with a typical uncertainty 
of 0.2 % at 30 and 44 GHz, and 0.1 % at 70 GHz, i.e., a factor of 
two better than the value achieved in 2013. 
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Fig. 1. Noise stationarity for a selection of two bolometers. The left panels show the total noise trends for each bolometer (dots). 
The solid line shows a running box average. The black dots are from the 2013 data release and the blue dots concern this release. The 
right panels show a histogram of the trends on the left. The box gives the width of the distribution at half maximum, as measured 
on the histogram, normalized to the mean noise level. The time response deconvolution has changed between the two data release 
and hence the absolute noise level is different. 
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Fig. 2. Fraction of discarded data per bolometer (squares with thick black line). The fraction of data discarded from glitch-flagging 
alone is shown with stars and the thin green line. The blue line with diamonds indicates the average fraction of discarded samples 
in valid rings. The two RTS bolometers (143_8 and 545_3) are not shown, since they are not used in the data processing. 
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The scanning beams^ used in the LFI pipeline (affecting 
calibration, effective beams, and beam window functions) are 
based on GRASP simulations, properly smeared to take into ac¬ 
count the satellite motion, and are similar to those presented in 
Planck Collaboration IV (2014). They come from a tuned optical 
model, and represent the most realistic fit to the available mea¬ 
surements of the LFI main beams. In Planck Collaboration IV 
(2014), calibration was performed assuming a pencil beam, the 
main beams were full-power main beams, and the resulting beam 
window functions were normalized to unity. For the 2015 re¬ 
lease, a different beam normalization has been used to prop¬ 
erly take into account the power entering the main beam (typ¬ 
ically about 99 % of the total power). Indeed, as described in 
Planck Collaboration V (2015), the current LFI calibration takes 
into account the full An beam (i.e., the main beam, as well as near 
and far sidelobes). Consequently, in the calculation of the win¬ 
dow function, the beams are not normalized to unity; instead, 
their normalization uses the value of the efficiency calculated 
taking into account the variation across the band of the optical 
response (coupling between feed horn pattern and telescope) and 
the radiometric response (band shape). 

Although the GRASP beams are computed as the far-held an¬ 
gular transmission function of a linearly polarized radiating ele¬ 
ment in the focal plane, the far-held pattern is in general not per¬ 
fectly linearly polarized, because there is a spurious component 
induced by the optical system, named “beam cross-polarization.” 
The Jupiter scans allowed us to measure only the total held, that 
is, the CO- and cross-polar components combined in quadrature. 
The adopted beam model has the added value of dehning the co- 
and cross-polar pattern separately, and it permits us to properly 
consider the beam cross-polarization in every step of the LFI 
pipeline. The GRASP model, together with the pointing informa¬ 
tion derived from the focal plane geometry reconstruction, gives 
the most advanced and precise noise-free representation of the 
LFI beams. 

The polarized main beam models were used to calculate the 
effective beams^, which take into account the specihc scanning 
strategy in order to include any smearing and orientation ef¬ 
fects on the beams themselves. Moreover, the sidelobes were 
used in the calibration pipeline to correctly evaluate the gains 
and to subtract Galactic straylight from the calibrated timelines 
(Planck Collaboration II 2015). 

To evaluate the beam window functions, we adopted two in¬ 
dependent approaches, both based on Monte Carlo simulations. 
In one case, we convolved a fiducial CMB signal with realistic 
scanning beams in harmonic space to generate the corresponding 
timelines and maps. In the other case, we convolved the fiducial 
CMB map with effective beams in pixel space with the FEBeCoP 
(Mitraetal. 2011) method. Using the first approach, we have 
also evaluated the contribution of the near and far sidelobes on 
the window functions. The impact of sidelobes on low mul¬ 
tipoles is about 0.1 % (for details see Planck Collaboration IV 
2015). 


^The term “scanning beam” refers to the angular response of a sin¬ 
gle detector to a compact source, including the optical beam and (for 
HFI) the effects of time domain filtering. In the case of HFI, a Fourier 
filter deconvolves the bolometer/electronics time response and lowpass- 
filters the data. In the case of LFI, the sampling tends to smear signal in 
the time domain. The term “effective beam” refers to a beam defined in 
the map domain, obtained by averaging the scanning beams pointing at 
a given pixel of the sky map taking into account the scanning strategy 
and the orientation of the beams themselves when they point along the 
direction to that pixel. See (Planck Collaboration IV 2014). 


The error budget was evaluated as in the 2013 release and 
it comes from two contributions: the propagation of the main 
beam uncertainties throughout the analysis; and the contribution 
of near and far sidelobes in the Monte Carlo simulation chain. 
The two error sources have different relevance, depending on 
the angular scale. Ignoring the near and far sidelobes is the dom¬ 
inant error at low multipoles, while the main beam uncertainties 
dominate the total error budget at ^ > 600. The total uncertain¬ 
ties in the effective beam window functions are 0.4 % and 1 % at 
30 and 44 GHz, respectively (at £ ^ 600), and 0.3 % at 70 GHz 
at ^ ^ 1000. 

5.2.2. HFI beams 

The HFI main beam measurement is described in detail 
in Planck Collaboration VII (2015) and is similar to that of 
Planck Collaboration VII (2014), although with several impor¬ 
tant changes. The HFI scanning beam model is a “Bspline” de¬ 
composition of the time-ordered data from planetary observa¬ 
tions. The domain of reconstruction of the main beam is ex¬ 
tended from a 40' square to a 100' square and is no longer 
apodized in order to preserve near sidelobe structure in the main 
beam model Planck Collaboration XXXI (2014), as well as to 
incorporate residual time-response effects into the beam model. 
A combination of Saturn and Jupiter data (rather than Mars 
data) is used for an improved signal-to-noise ratio, and a simple 
model of diffraction (consistent with physical optics predictions) 
is used to extend the beam model below the noise floor from 
planetary data. A second stage of cosmic ray glitch removal is 
added to reduce bias from unflagged cosmic ray hits. 

The effective beams and effective beam window functions 
are computed using the FEBeCoP and Quickbeam codes, as in 
Planck Collaboration VII (2014). While the scanning beam mea¬ 
surement produces a total intensity map only, effective beam 
window functions appropriate for both temperature and polar- 
ixed angular power spectra are produced by averaging the indi¬ 
vidual detector window functions weighed by temperature sen¬ 
sitivity and polarization sensitivity. Temperature-to-polarization 
leakage due to main beam mismatch is subdominant to noise in 
the polarization measurement, and is corrected as an additional 
nuisance parameter in the likelihood. 

The uncertainty in the beam measurement is derived from an 
ensemble of 100 Monte Carlo simulations of the planet obser¬ 
vations that include random realizations of detector noise, cos¬ 
mic ray hits, and pointing uncertainty propagated through the 
same pipeline as the data to simulated scanning beam products 
and simulated effective beam window functions. The error is ex¬ 
pressed in multipole space as a set of error eigenmodes, which 
capture the correlation structure of the errors. Additional con¬ 
sistency checks are performed to validate the error model, split¬ 
ting the planet data to construct Year 1 and Year 2 beams and to 
create Mars-based beams. With improved control of systemat- 
ics and higher signal-to-noise ratio, the uncertainties in the HFI 
beam window function have decreased by more than a factor of 
10 relative to the 2013 release. 

Several differences between the beams in 2013 and 2015 
may be listed. 

- Finer polar grid. Instead of the cartesian grid 40' on each 
side used previously, the beam maps were produced on both 
a cartesian grid of 200' on each side and 2" resolution, and 
a polar grid with a radius of 100' and a resolution of 2" in 
radius and 30' in azimuth. The latter grid has the advantage 
of not requiring any extra interpolation to compute the beam 
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spherical harmonic coefficients required by quickbeam, 
and therefore improves the accuracy of the resulting B{i). 

- Scanning beam elongation. To account for the elongation of 
the scanning beam induced by the time response deconvolu¬ 
tion residuals, the quickbeam computations are conducted 
with the birn for -6 < m < 6. We checked that the miss¬ 
ing terms account for less than 10“"^ of the effective at 
^ = 2000. Moreover, spotcheck comparisons with the effec¬ 
tive B(^) obtained by FEBeCoP show very good agreement. 

“ Finite size of Saturn . Even though its rings seem invisi¬ 
ble at Planck frequencies (and unlike Mars), Saturn has an 
angular size that must be accounted for in the beam win¬ 
dow function. The planet was assumed to be a top-hat disc 
of radius 9.5" at all HFI frequencies, whose window func¬ 
tion is well approximated by that of a 2D Gaussian profile of 
FWHM 1 F'185; the effective B{1) were therefore divided by 
that window function. 

“ Cut sky and pixel shape variability. The effective beam 
window functions do not include the (nominal) pixel win¬ 
dow function, which must be accounted for separately when 
analysing Planck maps. However, the shape and individual 
window function of the HEALPix Gorski et al. (2005) pix¬ 
els have large-scale variations around their nominal values 
across the sky. These variations impact the effective beam 
window functions applicable to Planck maps, in which the 
Galactic plane has been masked more or less conservatively, 
and are included in the effective B{f) that are provided. 

- Polarization and detector weights. Each 143, 217 and 
353 GHz frequency map is a combination of measurement by 
polarization-sensitive and polarization-insensitive detectors, 
each having a different optical response. As a consequence, 
at each of these frequencies, the Q and U maps will have a 
different beam window function than the I map. When cross- 
correlating the 143 and 217 GHz maps for example, the TT, 
EE, TE, and ET spectra will each have a different beam 
window function. 

- Polarization and beam mismatch. Since polarization mea¬ 
surements are differential by nature, any mismatch in 
the effective beams of the detectors involved will cou¬ 
ple with temperature anisotropies to create spurious polar¬ 
ization signals (Hu et al. 2003). In the likelihood pipeline 
(Planck Collaboration XI 2015) this additive leakage is mod¬ 
elled as a polynomial whose parameters are fit on the power 
spectra. 

“ Beam error model. See above. The improved S/N compared 
to 2013 leads to smaller uncertainties. At ^ = 1000 the un¬ 
certainties on B^ are 2.2 x 10“"^, 0.84 x 10“"^, and 0.81 x 10“"^ 
for 100, 143, and 217 GHz, respectively. At ^ = 2000, they 
are 11 X lO'^, 1.9 x lO'^, and 1.3 x lO'^. 

A reduced instrument model (RIMO) containing the effec¬ 
tive B(^) for temperature and polarization detector assemblies 
will be provided, for both auto- and cross-spectra. The RIMO 
will also contain the beam error eigenmodes and their covari¬ 
ance matrices. 

5.3. Focal plane geometry and pointing 

The focal plane geometry of LFI was determined independently 
for each Jupiter crossing (Planck Collaboration IV 2015), using 
the same procedure adopted in the 2013 release. The solutions 
for the seven crossings agree within 4" at 70 GHz (and 7" at 30 


and 44 GHz). The uncertainty in the determination of the main 
beam pointing directions evaluated from the single scans is about 
4" for the nominal scans, and 275 for the deep scans at 70 GHz 
(27" for the nominal scan and 19" for the deep scan, at 30 and 
44 GHz). Stacking the seven Jupiter transits, the uncertainty in 
the reconstructed main beam pointing directions becomes 076 
at 70 GHz and 2" at 30 and 44 GHz. With respect to the 2013 
release, we have found a difference in the main beam pointing 
directions of about 5" in the cross-scan direction and 076 in the 
in-scan direction. 

Throughout the extended mission, Planck continued to op¬ 
erate star camera STRl, with the redundant unit, STR2, only 
briefly swapped on for testing. No changes were made to the 
basic attitude reconstruction. We explored the possibility of up¬ 
dating the satellite dynamic model and using the fibre-optic gyro 
for additional high frequency attitude information. Neither pro¬ 
vided significant improvements to the pointing and were actually 
detrimental to overall pointing performance; however, they may 
become useful in future attempts to recover accurate pointing 
during the “unstable” periods. 

Attitude reconstruction delivers two quantities: the satellite 
body reference system attitude; and the angles between it and 
the principal axis reference system (so-called “tilt” or “wobble” 
angles). The tilt angles are needed to reconstruct the focal plane 
line-of-sight from the raw body reference frame attitude. For un¬ 
known reasons, the reconstructed tilt angles became irregular at 
the start of the LFI-only extension (cf. Fig. 3). Starting near day 
1000 after launch, the tilt angles began to indicate a drift that 
covered 1.'5 over about a month of operations. We found that 
the drift was not present in observed planet positions and we 
were therefore forced to abandon the reconstructed tilt angles 
and include the tilt correction into our ad hoc pointing correc¬ 
tion, PTCOR. 

We noticed that the most significant tilt angle corrections 
prior to the LFI extension tracked well the distance, Jsun, be¬ 
tween the Sun and Planck (see Fig. 3, bottom panel), so we 
decided to replace the spline fitting from 2013 with the use of 
the Solar distance as a fitting template. The fit was improved by 
adding a linear drift component and inserting breaks at events 
known to disturb the spacecraft thermal enviroment. In Fig. 4 
we show the co- and cross-scan pointing corrections and a selec¬ 
tion of planet position offsets after the correction was applied. 
The template-based pointing correction differs only marginally 
from the 2013 PTCOR, but an update was absolutely necessary 
to provide consistent, high fidelity pointing for the entire Planck 
mission. 

Finally we addressed the LFI radiometer electronics box 
assembly (REBA) interference that was observed in the 2013 
release by constructing, fitting, and subtracting another tem¬ 
plate from the REBA thermometry. This greatly reduced short 
timescale pointing errors prior to REBA thermal tuning on 
day 540 after launch. The REBA template removal reduced the 
pointing period timescale errors from 277 to 078 (in-scan) and 
179 (cross-scan). 

5.4. Calibration 

In this section we compare the relative photometric calibration of 
the all-sky CMB maps between LFI and HFI, as well as between 
Planck and WMAP. The two Planck instruments use different 
technologies and are subject to different foregrounds and sys¬ 
tematic effects. The Planck and WMAP measurements overlap 
in frequency range, but have independent spacecraft, telescopes, 
and scanning strategies. Consistency tests between these three 
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Fig. 3. Reconstructed tilt (wobble) angles between the satellite 
body frame and the principal axis frame. Vertical blue lines mark 
the ends of operation years and the dashed black line indicates 
day 540 after launch, when the thermal control on the LFI ra¬ 
diometer electronics box assembly (REBA) was adjusted. Top: 
First angle, , corresponds to a rotation about the satellite axis 
just 5° off the focal plane centre. Observed changes in only 
have a small impact on focal plane line-of-sight. Bottom: Second 
angle, is perpendicular to a plane defined by the nominal spin 
axis and the telescope line of sight. Rotation in 0^2 immediately 
impacts the opening angle and thus the cross-scan position of the 
focal plane. We also plot a scaled and translated version of the 
Solar distance that correlates well with if /2 until the reconstructed 
angles became compromised around day 1000 after launch. 

data sets are very demanding tests of the control of calibration, 
transfer functions, systematic effects, and foreground contami¬ 
nation. 

5.4.1. The orbital dipole 

In the 2013 data release, photometric calibration from 30 to 
353 GHz was based on the “Solar dipole”, that is, the dipole in¬ 
duced in the CMB by the motion of the Solar System barycentre 
with respect to the CMB. We used the value of the dipole mea¬ 
sured by WMAP5 (Hinshaw et al. 2009; Jarosik et al. 2011). 


Fig. 4 . Our ad hoc pointing correction, PTCOR, and a selection 
of observed planet position offsets after applying the correction. 
Top: Cross-scan pointing offset. This angle is directly affected 
by the second tilt angle, 02, in Fig. 3. Bottom: In-scan pointing 
offset. This angle corresponds to the spin phase and matches the 
third satellite tilt angle, 03 . Since 03 is poorly resolved by stan¬ 
dard attitude reconstruction, the in-scan pointing was already 
driven by PTCOR in the 2013 release. 

In the 2015 data release, photometric calibration of both LFI 
and HFI is based on the “orbital dipole,” i.e., the modulation of 
the Solar dipole induced by the orbital motion of the satellite 
around the Solar System barycentre. By using this primary cal¬ 
ibrator, we can derive for each Planck detector (or combination 
of detectors) an independent measurement of the Solar dipole, 
which is then used in the Planck calibration pipeline. The or¬ 
bital motion is known with exquisite accuracy, making the or¬ 
bital dipole potentially the most accurate calibration source in all 
of astrophysics, limited ultimately by the accuracy of the tem¬ 
perature of the CMB. The amplitude of this modulation, how¬ 
ever, is only about 250 yuK, varying with the details of the satel¬ 
lite motion, an order of magnitude smaller than the Solar dipole. 
Realizing its advantages as a fundamental calibration source re¬ 
quires low noise and good control of foregrounds, sidelobes, and 
large-angular-scale systematics. For the 2015 release, improve¬ 
ments in the control of systematic effects and foregrounds on 
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both LFI and HFI, including the availability of 2.5 and 4 orbital 
cycles for HFI and LFI, respectively (compared to 1.25 cycles 
in the 2013 release), have allowed accurate calibration of both 
instruments on the orbital dipole, summarized in the following 
subsections and described in detail in Planck Collaboration II 
(2015) and Planck Collaboration VIII (2015). The dipole com¬ 
ponent of the CMB and the frequency-independent part of the 
quadrupole (induced by the Solar dipole) are removed from both 
the LFI and HFI data; however, higher order effects of the Solar 
dipole (see Planck Collaboration XXVII 2014) are left in the 
data, which is matched by what is contained in the simulations 
Planck Collaboration XII (2015). 

With the 2015 data calibrated on the orbital dipole, Planck 
has made independent measurements of the Solar dipole 
(Table 1), which can be compared to the WMAP5 measurement 
(Hinshaw et al. 2009). Amplitudes agree within 0.28%; direc¬ 
tions agree to better than 2'. Although the difference in am¬ 
plitude between the Planck and the WMAP5 measurements of 
the Solar dipole is small and within uncertainties, it had non- 
negligible consequences in 2013. WMAP was calibrated on the 
orbital dipole, so errors in its Solar dipole measurement did 
not contribute to its overall calibration errors. Planck in 2013, 
however, was calibrated on the WMAP5 Solar dipole, which 
is 0.28% lower than the orbital-dipole-calibrated 2015 Planck 
measurement. Calibrating LFI and HFI against WMAP5 in the 
2013 results, therefore, resulted in 2013 gains that were 0.28 % 
too low for both LFI and HFI. This factor is included in Tables 2 
and 3. 

5.4.2. Instrument level calibration 

LFI —There were four significant changes related to LFI cali¬ 
bration between the 2013 and the 2015 results. First (as antici¬ 
pated in the 2013 LFI calibration paper, Planck Collaboration V 
2014), the convolution of the beam with the overall dipole (Solar 
and orbital dipoles, including their induced frequency indepen¬ 
dent quadrupoles) is performed with the full An beam rather 
than a pencil beam. This dipole model is used to extract the 
gain calibration parameter. Because the details of the beam pat¬ 
tern are unique for each detector even within the same fre¬ 
quency channel, the reference signal used for the calibration is 
different for each of the 22 LFI radiometers. This change im¬ 
proves the results of null tests and the quality of the polarization 
maps. When taking into account the proper window function 
(Planck Collaboration IV 2015), the new convolution scheme 
leads to a shifts of -rO.32, -r0.03, and -rO.30% in gain calibra¬ 
tion at 30, 44, and 70 GHz, respectively (see Table 2). Second, 
a new destriping code. Da Capo (Planck Collaboration V 2015), 
is used; this supersedes the combination of a dipole fitting rou¬ 
tine and the Mademoiselle code used in the 2013 data release 
and offers improved handing of 1 // noise and residual Galactic 
signals. Third, Galactic contamination entering via sidelobes is 
subtracted from the timelines after calibration. Finally, a new 
smoothing algorithm is applied to the calibration parameters. It 
adapts the length of the smoothing window depending on a num¬ 
ber of parameters, including the peak-to-peak amplitude of the 
dipole seen within each ring and sudden temperature changes in 
the instrument. These changes improve the results of null tests, 
and also lead to overall shifts in gain calibration a few per mill, 
depending on frequency channel. The values reported in the third 
column of Table 2 are approximate estimates from the combi¬ 
nation of improved destriping. Galactic contamination removal, 
and smoothing. They are calculated under the simplifying as¬ 
sumption that these effects are completely independent of the 


beam convolution and can therefore be combined linearly with 
the latter (for more details see Planck Collaboration V (2015)). 

In total, these four improvements give an overall increase in 
gain calibration for LFI of -r0.17, -rO.36, and -rO.54 % at 30, 44, 
and 70 GHz, respectively. Adding the 0.28 % error introduced by 
the WMAP Solar dipole in 2013 (discussed in Sect. 5.4.1), for 
the three LFI frequency channels we find overall shifts of about 
0.5, 0.6 and 0.8 % in gain calibration with respect to our LFI 
2013 analysis (see Table 2). 

As shown in Planck Collaboration V (2015), relative calibra¬ 
tion between LFI radiometer pairs is consistent within their sta¬ 
tistical uncertainties. At 70 GHz, using the deviations of the cal¬ 
ibration of single channels, we estimate that the relative calibra¬ 
tion error is about 0.10 %. 

HFI —There were three significant changes related to HFI cal¬ 
ibration between the 2013 and the 2015 results: improved de¬ 
termination and handling of near and far sidelobes; improved 
ADC non-linearity correction; and improved handling of very 
long time constants. The most significant changes arise from the 
introduction of near sidelobes of the main beam (0?5-5°) that 
were not detected in observations of Mars, and from the intro¬ 
duction of very long time constants. We consider these in turn. 

Observations of Jupiter were not used in 2013 results, be¬ 
cause its signal is so strong that it saturates some stages of the 
readout electronics. The overall transfer function for each de¬ 
tector is corrected through the deconvolution of a time transfer 
function, leaving a compact effective beam used together with 
the maps in the science analysis. In the subsequent “consis¬ 
tency paper” Planck Collaboration XXXI (2014), it was found 
that lower-noise hybrid beams built by using Mars, Saturn, and 
Jupiter observations reveal near sidelobes leading to significant 
corrections of 0.1 to 0.3 %. Far sidelobes give a very small cal¬ 
ibration correction that is almost constant for ^ > 3. The zodi¬ 
acal contribution was removed in the timelines, as it does not 
project properly on the sky. It gives an even smaller and negligi¬ 
ble correction, except in the submillimetre (hereafter “submm”) 
channels at 545 and 857 GHz. 

The most significant change results from the recognition of 
the existence of very long time constants (VLTC) and their in¬ 
clusion in the analysis. VLTCs introduce a significant shift in the 
apparent position of the dominant anisotropy in the CMB, the 
Solar dipole, away from its true position. This in effect creates 
a leakage of the Solar dipole into the orbital dipole. This distur¬ 
bance is the reason why calibration on the orbital dipole did not 
work as expected from simulations, and why calibration in 2013 
was instead based on the WMAP5 Solar dipole. As discussed 
in Sect. 5.4.1, the WMAP5 Solar dipole was underestimated by 
0.28 % when compared with the Planck best-measured ampli¬ 
tude, leading to an under-calibration of 0.28 % in the Planck 
2013 maps. With VLTCs included in the analysis, calibration 
on the orbital dipole worked as expected, and gave more accu¬ 
rate results, while at the same time eliminating the need for the 
WMAP5 Solar dipole and removing the 0.28 % error that it in¬ 
troduced in 2013. 

These HFI calibration changes are reported in Table 3. 
Together, they give an average shift of gain calibration of typ¬ 
ically 1 % (Planck Collaboration VIII 2015) for the three CMB 
channels, accounting for the previously unexplained difference 
in calibration on the first acoustic peak observed between HFI 
and WMAP. 

The relative calibration between detectors operating at the 
same frequency is within 0.05 % for 100 and 143 GHz, 0.1 % 
at 217 GHz, and 0.4% at 353 GHz Planck Collaboration VIII 
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Table 1. LFI, HFI, and WMAP measurements of the Solar dipole. 


Galactic coordinates 


Amplitude I b 

Experiment [//Kcmb] [deg] [deg] 

LFP. 3365.5 +3.0 264.01 +0.05 48.26 +0.02 

HFP. 3364.29 + 1.1 263.914 + 0.013 48.265 + 0.002 

Planck 2015 nominaP. 3364.5 + 2.0^ 264.00 +0.03 48.24 +0.02 

WMAP^ . 3355^ + 8 263.99 +0.14 48.26 +0.03 


^ The “nominal” Planck dipole was chosen as a plausible combination of the LFI and HFI measurements early in the analysis, to carry out 
subtraction of the dipole from the frequency maps (see Sect. 5.4.3). The current best determination of the dipole comes from an average of 100 
and 143 GHz results (Planck Collaboration VIII 2015). 

^ Uncertainties include an estimate of systematic errors. 

^ Hinshaw et al. (2009). 

^ See Sect. 5.4.1 for the effect of this amplitude on Planck calibration in 2013. 


Table 2. LFI calibration changes at map level, 2013 ^ 2015. 


Beam solid Pipeline Orbital 


Frequency angle improvements^ Dipole*^ Total 

[GHz] [%] [%] [%] [%] 

30 . 0.32 -0.15 0.28 0.45 

44 . 0.03 +0.33 0.28 0.64 

70. +0.30 +0.24 0.28 0.82 


^ This term includes the combined effect of the new destriping code, 
subtraction of Galactic contamination from timelines, new smoothing 
algorithm. It has been calculated under the hypothesis that it is fully 
independent of the beam convolution. 

^ Change from not being dependent on the amplitude error of the 
WMAP9 Solar dipole (Sect. 5.4.1). 

Table 3. HFI calibration changes at map level, 2013 ^ 2015. 


SiDELOBES Orbital dipole 


Frequency Near Far Dipole^ VLTC Total 

[GHz] [%] [%] [%] [%] [%] 


100 . 0.2 0.087 0.28 0.49 1.06 

143 . 0.2 0.046 0.28 0.47 1.00 

217 . 0.2 0.043 0.28 0.66 1.17 

353 . 0.275 0.006 0.28 1.5 2.06 


^ Change from not being dependent on the amplitude error of the 
WMAP9 Solar dipole (Sect. 5.4.1). 

(2015). These levels for the CMB channels are within a fac¬ 
tor of 3 of the accuracy floor set by noise in low T’ polarization 
(Tristram et al. 201 1). 

The 545 and 857 GHz channels are calibrated separately us¬ 
ing models of planetary atmospheric emission. As in 2013, we 
used both Neptune and Uranus. The main difference comes from 
better handling of the systematic errors affecting the planet flux 
density measurements. Analysis is now performed on the time¬ 
lines, using aperture photometry, and taking into account the in¬ 
homogeneous spatial distribution of the samples. On frequency 
maps, we estimate statistical errors on absolute calibration of 
1.1 % and 1.4 % at 545 and 857 GHz, respectively, to which we 
add the 5 % systematic uncertainty arising from the planet mod¬ 
els. Errors on absolute calibration are therefore 6.1 and 6.4% 


at 545 and 857 GHz, respectively. Since the reported relative un¬ 
certainty of the models is of the order of 2 %, we And the relative 
calibration between the two HFI high-end frequencies to be bet¬ 
ter than 3 %. Relative calibration based on diffuse foreground 
component separation gives consistent numbers (see table 6 of 
Planck Collaboration X 2015). Compared to 2013, calibration 
factors changed by 1.9 and 4.1 %, at 545 and 857 GHz, respec¬ 
tively. Combined with other pipeline changes (such as the ADC 
corrections), the brightness of the released 2015 frequency maps 
has decreased by 1.8 and 3.3 % compared to 2013. 


5.4.3. Relative calibration and consistency 

The relative calibration of LFI, HFI, and WMAP can be assessed 
on several angular scales. At 7’ = 1, we can compare the ampli¬ 
tude and direction of the Solar dipole, as measured in the fre¬ 
quency maps of the three instruments. On smaller scales, we can 
compare the amplitude of the CMB fluctuations measured fre¬ 
quency by frequency by the three instruments, during and after 
component separation. 

Comparison of independent measurements of the Solar 
dipole—Table 1 gives the LFI and HFI measurements of the 
Solar dipole, showing agreement well within the uncertainties. 
The amplitudes agree within 1.5 juK (0.05 %), and the directions 
agree within 3'. Table 1 also gives the “nominal” Planck dipole 
that has been subtracted from the Planck frequency maps in the 
2015 release. This is a plausible combination of the LFI and HFI 
values, which satisfied the need for a dipole that could be sub¬ 
tracted uniformly across all Planck frequencies early in the data 
processing, before the Anal systematic uncertainties in the dipole 
measurements were available and a rigorous combination could 
be determined. See Planck Collaboration VIII (2015) Sect. 5.1 
for additional measurements. 

Nearly independent determinations of the Solar dipole can 
be extracted from individual frequency maps using component- 
separation methods relying on templates from low and high fre¬ 
quencies where foregrounds dominate (Planck Collaboration V 
2015; Planck Collaboration VIII 2015). The amplitude and di¬ 
rection of these Solar dipole measurements can be compared 
with each other and with the statistical errors. This leads to rela¬ 
tive gain calibration factors for the 7’ = 1 mode of the maps ex¬ 
pressed in Kcmb, as shown for frequencies from 70 to 545 GHz 
in Table 4. For components of the signal with spectral distribu- 
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tion different from the CMB, a colour correction is needed to 
take into account the broad bands of these experiments. 

Comparison of the residuals of the Solar dipole left in the 
CMB maps after removal of the best common estimate— 

Another measurement of relative calibration is given by the 
residuals of the Solar dipole left in CMB maps after remov¬ 
ing the best common estimate, the nominal Planck dipole. The 
residual dipole comes from two terms as illustrated in Fig. 5, 
one associated with the error in direction, with an axis nearly 
orthogonal to the Solar dipole, and one associated with the error 
in amplitude aligned with the Solar dipole. Using the 857 GHz 
map as a dust template (extrapolated with optimized coefficients 
derived per patch of sky), we find residual dipoles dominated by 
errors orthogonal to the direction of the Solar dipole at 100 and 
143 GHz, and residuals associated with calibration errors for the 
other frequencies. The relative residual amplitudes are given in 
Table 4. This shows that a minimization of the dipole residuals 
can and will be introduced in the HFI calibration pipeline for the 
final 2016 release. 



0.01 0.10 1.00 10.00 
6G [%] 


Fig. 5. Angle difference Ofres-rm between the removed Solar 
dipole and the residual dipole for given errors on the dipole di¬ 
rection (i.e., the angle difference between the removed dipole 
and the true Solar dipole, (Tsoi-rm) and on the calibration {6G = 
1 - ^rmMsob expressed in percent). 


Comparison of CMB anisotropies frequency by fre¬ 
quency during and after component separation—Table 4 
also shows the relative calibration between frequencies and 
detectors determined by SMICA (Planck Collaboration XV 
2014; Planck Collaboration IX 2015) and Commander 
(Planck Collaboration IX 2015; Planck Collaboration X 2015), 
two of the map-based diffuse component separation codes used 
by Planck. The calculation is over different multipole ranges 
for the two methods, so variation between the two could reflect 
uncertainties in transfer functions; moreover. Commander uses 
different constraints in order to deal with the complexities and 
extra degrees of freedom involved in fitting foregrounds indi¬ 
vidually (see Planck Collaboration X 2015 for details), so we 
do not expect identical results with the two codes. Nevertheless, 
the agreement is excellent, at the 0.2 % level between the first 


acoustic peak, intermediate i, and dipole residuals, and the 
intercalibration offsets between frequencies are within 0.3 % of 
zero from 30 GHz to 217 GHz. 

The following points highlight the remarkable internal con¬ 
sistency of Planck calibration. 

- The small Solar dipole residuals measured for the 100 and 
143 GHz channels (< 4yuK) are close to 90°away from the 
adopted Planck Solar dipole, reflecting in both cases a small 
2.'8 shift in the measured direction of the dipole compared 
to the adopted dipole, but amplitudes (hence calibrations) 
within 0.1 % of the adopted (“mean”) value. The Commander 
and SMICA inter-comparisons below and on the first acoustic 
peak give a calibration difference between 100 and 143 GHz 
of < 0.09 %, confirming the very high calibration accuracy 
of these two channels. 

- The amplitude of the Solar dipole measured by the 70 GHz 
channel shows a difference of 1 yuK (0.03 %) with respect to 
the best HFI Solar dipole amplitude. 

- The 217, 353, and 545 GHz channels show dipole residuals 
aligned with the Solar dipole, which thus measure directly 
calibration errors with respect to 143 GHz of 0.2, 0.53, and 
1.25%. 

- The SMICA first peak intercalibration of 217 and 353 GHz 
with respect to 143 GHz, taken again as reference, shows 
similar intercalibration to the dipole residuals with differ¬ 
ences 0.08 and 0.20 %. In fact. Table 4 suggests that we can 
now achieve significantly better intercalibration of all CMB 
channels from 70 to 353 GHz. 

- Comparison of the Solar dipole and first acoustic peak inter¬ 
calibration factors for the 545 GHz channel gives a difference 
of only 0.16 %. This shows that the 545 GHz channel could 
be calibrated using the first acoustic peak of the CMB instead 
of planets. Use of the planet model could then be limited to 
intercalibration between 545 and 857 GHz. The roughly 1 % 
agreement between the planets and CMB absolute calibra¬ 
tions also show that the current uncertainties on the absolute 
calibration of the high-frequency channels, dominated by the 
roughly 5 % error on the models, are probably overestimated. 

- The intercalibration factors derived from Commander in 
all frequency bands from 70 GHz to 217 GHz are less 
than 0.1 %. Considering all Planck bands from 30 GHz to 
353 GHz, they are within 0.5 %. 

This comparison can also be made at the power spectrum 
level. The left-hand panel in Fig. 6 compares the 70, 100, and 
143 GHz channels of LFI and HFI in the multipole range of the 
first acoustic peak, 50 < ^ < 500, uncorrected for foregrounds, 
over 60 % of the sky. The low values a ^ = 50 show the effect 
of unremoved diffuse foregrounds at 143 GHz, and the rise of 
the 70/143 ratio is at least partly driven by unremoved discrete 
foregrounds; the uncertainties are larger at 70 GHz as well. In 
the middle region, the agreement is very good, at a level of a few 
tenths of a percent. This result is a direct test that all systematic 
effects in calibration have been corrected on both instruments to 
better than this value. 

The right-hand panel of Fig. 6 shows the ratios of Planck 
TT spectra at 70 and 100 GHz to those of WMAP in 
the V and W bands, as calculated for Planck 2013 data 
(Planck Collaboration XXXI2014) and for the 2015 data. While 
the scatter is significantly larger than that in the left-hand panel, 
due to the higher noise in WMAP, the agreement is very good, 
and within the statistical errors. We can now say that within 
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Table 4. Intercalibration factors by frequency between LFI, HFI, and WMAR 


CMB Anisotropy [%] 


Frequency [GHz] (Detector) 

Solar dipole [%] 

1 

Commander 
25<i< 100 

SMICA 

50 < ^ < 500 

30. 


-0.3" +0.1 


44. 


0.3" +0.1 


70. 

0.04" 

0.0" +0.1 

0.21"+ 0.06 

100. 

0.03 

0.09 + 0.02 

0.03 + 0.02 

143. 

0*^ 

0; -0.1 +0.F 

0^ 

217. 

0.20 

0; 0.02 + 0.03^ 

0.28 + 0.02 

353 . 

0.53 

0.5 +0.1 

0.73 + 0.11 

545 . 

1.25 

-1.0^ 

1.09 + 1.5 

WMAP (23) . 


0^ 


WMAP (33) . 


0.1 +0.1 


WMAP 41 (Q). 


0.1 +0.1 


WMAP 61 (V). 


0.2 +0.1 


WMAP 94 (W). 

-0.26 

0.2 +0.1 

0.28 + 0.15 


^ LFI map rescaling factors that are incorporated in the beam transfer functions, as described in Planck Collaboration II (2015), have been applied. 

^ Reference frequency; no intercalibration calculated. 

^ For Commander at 143 GHz, detector set ‘dsl’ was used as a reference (intercalibration factor = 0). The mean intercalibration factor for detectors 
ds2+5+6+7 was -0.1 +0.1. Similarly, at 217 GHz detector ‘1’ was used as a reference (intercalibration factor = 0), and the mean intercalibration 
factor for detectors 2+3+4 was 0.02 + 0.03. See Table 6 in Planck Collaboration X (2015) for details. 

^ For Commander, the effective recalibration of the 545 GHz channel measured in units of yuKcmb is the product of a multiplicative calibration 
factor and a unit conversion correction due to revised bandpass estimation. See Sect. 5.3 in Planck Collaboration X (2015) for details. 



Fig. 6. Ratios of power spectra over the region of the first acoustic peak, uncorrected for foregrounds (which vary over the three 
frequencies), over 60 % of the sky. The uncertainties are the errors on the mean within each M = 40 bin of the ratios computed I 
by L Left: Ratios of 70 and 100 GHz TT spectra to 143 GHz. The low values at 7’ = 50 are due to diffuse foregrounds at 143 GHz. 
The rise to higher multipoles in the 70/143 ratio is due to discrete foregrounds. Right: Ratio of TT spectra of Planck 70 and 
100 GHz to WMAP V and W bands, as calculated for Planck 2013 data (Planck Collaboration XXXI 2014) and for the 2015 data. 
The near-overlap of frequencies between the Planck and WMAP bands means that foregrounds have no appreciable effect on the 
ratios. The effect of the calibration changes in Planck between 2013 and 2015 that are discussed in this paper is clear. There is now 
excellent agreement within statistical errors between Planck and WMAP in the region of the spectrum where both have high S/N. 
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the uncertainties, LFI, HFI, and WMAP agree, and the differ¬ 
ence seen in the 2013 data (Planck Collaboration XXXI 2014) 
is gone. 

5.4.4. Summary of calibration 

The Planck 70 and 100 GHz channels belong to instruments 
based on different technologies, with different systematic effects, 
and close to the minimum of the diffuse foregrounds. They thus 
provide a very good test of the consistency of calibration and 
transfer functions. The internal consistency between LFI and 
HFI is remarkable. Figure 6 represents a stringent test of cal¬ 
ibration, systematic effects, beams, and transfer functions, and 
demonstrates overall consistency at a level of a few parts per 
thousand between independent instruments and spacecraft. 

The Planck CMB-channels from 70 to 217 GHz show cal¬ 
ibration difference below 0.3 %, measured from both residual 
dipoles and first acoustic peak. Using a Solar dipole reference 
established on the 100 and 143 GHz channels, it is likely that 
all detectors could be inter-calibrated to 0.05 % in subsequent 
data processing versions. The agreement of the measured cali¬ 
bration factors from dipole residuals {i = \) and first acoustic 
peak {i = 200) shows that the transfer functions are controlled 
to better than 0.2 % in this multipole range. Corrections for sys¬ 
tematic effects in HFI cover a dynamic range from detector to 
detector larger than 2 at 100 and 143 GHz, but have reduced the 
calibration errors by an order of magnitude. This suggests that 
the corrections lead now to an absolute photometric calibration 
accuracy on the orbital dipole (limited only by systematics and 
noise) of 0.1 %. 

As in other instances in the Planck data processing, when 
very small systematic effects are detected and measured in a 
posteriori characterization, their removal from the data is com¬ 
plicated. Their determinations are often degenerate, and com¬ 
plete reprocessing is necessary. The calibration improvement 
demonstrated by the minimization of the dipole residuals us¬ 
ing the 857 GHz dust template will be introduced in a self- 
consistent way in the HFI calibration pipeline and overall pro¬ 
cessing for the final 2016 release. Furthermore, the use of the 
Solar dipole parameters from the best Planck CMB channels 
(100 and 143 GHz) will be introduced in the processing of the 
channels more affected by foregrounds and noise. The LFI cali¬ 
bration accuracy is now close to noise-limited, but improvements 
will be made in 2015 according to a complete simulation plan to 
improve our understanding of calibration and systematics affect¬ 
ing low multipoles, particularly for polarization analysis. 


6. Timelines 

For the first time, the 2015 Planck release includes time series of 
the observations acquired by individual detectors in LFI and HFI 
(see Planck Collaboration II 2015; Planck Collaboration VII 
2015, for details). These timelines will be of use for those wish¬ 
ing to construct maps using specific time periods or mapmaking 
algorithms. 

The delivered timelines have been cleaned of all major in¬ 
strumental systematic effects. For HFI timelines this means that 
the raw timelines are ADC-corrected, demodulated, despiked, 
corrected for rare baseline jumps, and a dark template has been 
removed; they are converted to absorbed power units, and the 
time transfer function has been deconvolved. For LFI timelines 
this means that the raw timelines are ADC-corrected, despiked, 
and demodulated; furthermore, the raw diode outputs (two per 


receiver) are combined and gain regularization is applied before 
calibration. 

The timelines are calibrated to astrophysical units and cor¬ 
rected for a zero-point level (determined at map level). The Solar 
and orbital dipole signals have been removed. In addition, for 
LFI, an estimation of Galactic straylight has been removed. 

The timelines still contain the low-frequency noise that is 
later removed by destriping at the mapmaking stage. However, 
a set of offsets are provided (determined during mapmaking), 
which can be used to convert the calibrated timelines to maps 
without destriping. For HFI a single offset per ring is determined; 
for LFI the offsets are computed every 0.246, 0.988, and 1.000 s 
for the 30, 44, and 70 GHz channels, respectively. For HFI, the 
offsets are determined during mapmaking using the full mission 
data set and all valid detectors per channel, and they are then 
applied to all the maps produced, i.e., using any fraction of the 
mission (year, survey) or any subset of detectors (single detector, 
detector set). For LFI the offsets are similarly determined using 
the full mission and all valid detectors per channel. These offsets 
have been used to produce the full-mission LFI maps; however, 
for shorter period maps, different offsets are used, which opti¬ 
mize noise cross-correlation effects, and these are not delivered. 

The timelines are accompanied by flags that determine which 
data have been used for mapmaking, as well as pointing time¬ 
lines, which are sampled at the same frequency as the data them¬ 
selves. 


7, Frequency Maps 

In Fig. 7 and Fig. 8 we show the Planck 2015 maps. Note that 
Planck uses HEALPix (Gorski et al. 2005) as its basic representa¬ 
tion scheme for maps, with resolution labelled by the Aside value. 

7.L Mapmaking 
7.1.1. LFI 

Mapmaking takes as its input the calibrated timelines, from 
which the cosmological and orbital dipole signals have been re¬ 
moved. An estimate of Galactic straylight is subtracted from the 
timelines prior to mapmaking, since this is difficult to correct 
for at map level. As for the 2013 release, the LFI maps are pro¬ 
duced using the Madam destriping code (Keihanen et al. 2010), 
enhanced with a noise prior, which enables accurate removal of 
correlated 1 // noise, while simultaneously minimizing system¬ 
atic errors by judicious use of masks. The production of maps 
and covariance matrices is validated using the FFP8 simulations. 
The output of the code consists of sky maps of temperature and 
Stokes Q and U polarization, and a statistical description of 
residual noise in the maps in the form of pixel-pixel noise co- 
variance matrices. These matrices are produced at Aside = 64. 
In addition to full-mission maps at both high and low resolution 
(^^side = 16), many other types of maps are produced, including 
those from single horns, single radiometers, single surveys, odd 
and even surveys, single years, and halves of the mission. The 
LFI maps are not corrected for beam shape, so that point sources 
in the map have the shape of the effective beam at that location. 
Zero-levels are estimated by fitting a cosecant-law Galactic lat¬ 
itude model to the CMB-subtracted maps, and subtracting this 
from the maps. The polarization maps must be corrected for 
bandpass leakage through multiplication with leakage template 
maps, which are estimated via a process similar to component 
separation. 
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Fig. 7. The nine Planck frequency maps show the broad frequency response of the individual channels. The color scale (identical to the one used in 2013), based on inversion of 
the function y = 10^ - 10"^, is tailored to show the full dynamic range of the maps. 
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Fig. 8. The seven Planck polarization maps between 30 and 353 GHz, shown in Stokes Q and U, as well as in total polarized intensity 
{P). The LFI maps are not bandpass-corrected; the HFI maps are. The color scale uses the same function as in Fig. 7, but the range 
limits have been adjusted. 


16 




















































Planck Collaboration: The Planck mission 


A summary of the characteristics of the LFI maps is pre¬ 
sented in Table 5. 

7.1.2. HFI 

As for the Planck 2013 release, the measurements in each 
HEALPix pixel visited during a stable pointing period (i.e., 
“ring”) are averaged for each detector, keeping track of the 
bolometer orientation on the sky. The calibration and mapmak¬ 
ing operations use this intermediate product as an input. For each 
detector, the TOIs are only modified by a single offset value 
per ring, determined using the destriping method described in 
Tristram et al. (2011). The offsets are computed simultaneously 
for all bolometers at a given frequency, using the full mission 
data. For a given bolometer, the same offset per ring is applied 
whatever the map (e.g., full-mission, half-mission, detector- 
set maps; but for half-ring maps, see Planck Collaboration VII 
2015). Each data sample is calibrated in ATcmb for the 100, 
143, 217, and 353 GHz channels, and MJy sr“^ (assuming vly = 
constant) for the 545 and 857 GHz channels, using the calibra¬ 
tion scheme presented in Sect. 5.4.2. Contrary to the 2013 re¬ 
lease, the bolometer gains are assumed to be constant through¬ 
out the mission. The final mapmaking is a simple projection of 
each unflagged sample to the nearest grid pixel. For polarization 
data, when several detectors are solved simultaneously, the po¬ 
larization mapmaking equation is inverted on a per-pixel basis 
(Planck Collaboration VII 2015). 

The products of the HFI mapmaking pipelines are pixelized 
maps of /, 2, and U, together with their covariances. Map res¬ 
olution is Aside = 2048, and the pixel size is 1.'07. The basic 
characteristics of the maps are given in Table 6. For details, see 
(Planck Collaboration VII 2015). 

Maps are cleaned for the zodiacal light component, which 
varies in time, based on templates fitted on the survey difference 
maps (see Planck Collaboration XIV 2014). These templates are 
systematically subtracted prior to mapmaking. The Planck total 
dipole (Solar and orbital) is computed and also subtracted from 
the data. Contrary to 2013, the far sidelobes (FSL) are not re¬ 
moved from the maps. 

The 2015 HFI maps delivered via the PL A have had zodiacal 
light removed, include the CIB and the zero level of the tempera¬ 
ture maps has been adjusted for Galactic emission. However, the 
zero level adjustment was based on maps which contained zodi¬ 
acal light, and therefore the released maps require an additional 
frequency-dependent correction which has to be applied man¬ 
ually. For work requiring all astrophysical sources of emission 
to be present in the maps, the corrections provided under Note 
“e2” of Table 6 must be added to the maps. For work requiring 
Galactic emission only, the “e2” corrections should be added to 
the maps, and the CIB levels provided under Note “e” should be 
removed. 

8. CMB Products 

8.1. CMB maps 

As for the Planck 2013 release, we use four different meth¬ 
ods to separate the Planck 2015 frequency maps into physi¬ 
cal components (Planck Collaboration IX 2015). The four meth¬ 
ods are: SMICA (independent component analysis of power 
spectra, Delabrouille et al. 2003; Cardoso etal. 2008); NILC 
(needlet-based internal linear combination, Delabrouille et al. 
2009); Commander (pixel-based parameter and template fitting 
with Gibbs sampling, Eriksen et al. 2006; Eriksen et al. 2008); 


and SEVEN (template fitting, Eernandez-Cobos et al. 2012). The 
methods used are conceptually the same as in 2013, but we ap¬ 
ply them now independently to the temperature and polarization 
maps. Similarly to what was done in 2013, simulations (in this 
case EEP8, Planck Collaboration XII 2015) are used to test the 
methods and estimate uncertainties in the recovery of compo¬ 
nents. 

All four methods produce CMB maps in Stokes /, 2, and 
U. In addition. Commander and SMICA also separate diffuse as¬ 
trophysical “foregrounds” characterized by their different spec¬ 
tral signatures. Commander does so by fitting physical models of 
these foregrounds and the CMB to the sky, whereas SMICA ex¬ 
tracts a fixed set of independent components representing CMB, 
foregrounds and noise; typically, SMICA assumes two “fore¬ 
grounds” are present at low and high frequencies, respectively. 
An important change in the implementation of Commander in 
2015 is the number of input maps used: firstly, the number of 
Planck maps is expanded to use detector-level maps rather than 
maps which combine all detectors at each frequency; secondly, 
the inputs include a map of 408 MHz emission, and the 9-year 
WMAP maps. The significant increase in the number of in¬ 
put maps allows Commander to: (a) control much better factors 
such as relative calibration and frequency response of individual 
channels; and (b) extract a larger number of foreground temper¬ 
ature components, now matching those that are expected to be 
present in the sky. 

The 2015 Planck CMB temperature maps produced by 
all four methods (see an example in Eig. 9) are signifi¬ 
cantly more sensitive than those produced in 2013 (by a fac¬ 
tor of 1.3). They are used mainly for non-Gaussianity analy¬ 
sis (Planck Collaboration XVII2015; Planck Collaboration XVI 
2015) and for the extraction of lensing deflection maps 
(Planck Collaboration XV 2015). Eor these analyses, all four 
methods are considered to give equivalently robust results, and 
the dispersion between the four gives a reasonable estimate of 
the uncertainty of the CMB recovery. Although the statistical 
properties of these maps give good results when used to fit cos¬ 
mological models, the best Planck 2015 cosmological parame¬ 
ters are derived from a likelihood code which allows much more 
detailed tuning of the contribution of individual frequencies 
and i-hy-i removal of foregrounds (Planck Collaboration XIII 
2015). A low-resolution version of the Commander map is also 
used in the pixel-based low-^ likelihood used to extract our best- 
fit 2015 cosmology (Planck Collaboration XI 2015). 

In polarization, the CMB maps resulting from the 2015 
Planck component separation methods represent a dramatic 
advance in terms of coverage, angular resolution, and sen¬ 
sitivity. Nonetheless, they suffer from a significantly high 
level of anomalous features at large angular scales, aris¬ 
ing from corresponding systematic effects in the input fre¬ 
quency maps between 100 and 217 GHz. The characteriza¬ 
tion of these systematic effects is still ongoing, and it is cur¬ 
rently suspected that low-level spurious signals are also present 
at intermediate angular scales (Planck Collaboration VII 2015; 
Planck Collaboration VIII 2015). Eor this reason, the CMB po¬ 
larization maps presented here have had their large angular 
scales ^ 30) filtered out. Eiltered Stokes Q and U maps re¬ 

sulting from one of the methods are shown as an example in 
Eig. 10. They are used only for a very limited number of cosmo¬ 
logical analyses, which have been shown to be immune to the 
undesired features still present: estimation of primordial non- 
Gaussianity levels (Planck Collaboration XVII 2015); stacking 
analysis (Planck Collaboration XVI 2015); estimation of pri¬ 
mordial magnetic field levels (Planck Collaboration XIX 2015); 
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Table 5. Main characteristics of LFI full mission maps. 


Characteristic 

30 GHz 

Frequency band 

44 GHz 

70 GHz 

Centre frequency [GHz]. 

28.4 

44.1 

70.4 

Effective beam FWHM^ [arcmin] . 

32.29 

27.00 

13.21 

Effective beam ellipticity^. 

1.32 

1.04 

1.22 

Temperature noise (1°)^ [^uKcmb]. 

2.5 

2.7 

3.5 

Polarization noise (1°)*^ [^uKcmb]. 

3.5 

4.0 

5.0 

Overall calibration uncertainty^ [%] . 

0.35 

0.26 

0.20 

Systematic effects uncertainty in Stokes [^uKcmb]. 

0.19 

0.39 

0.40 

Systematic effects uncertainty in Stokes [^uKcmb] • • • • 

0.20 

0.23 

0.45 

Systematic effects uncertainty in Stokes [^uKcmb] • • • • 

0.40 

0.45 

0.44 


^ Calculated from the main beam solid angle of the effective beam, flgff = mean(Q). These values are used in the source extraction pipeline 
(Planck Collaboration XXVI 2015). 

^ Noise rms computed after smoothing to 1°. 

^ Sum of the error determined from the absolute and relative calibration, see Planck Collaboration IV (2015). 

^ Estimated rms values over the full sky and after full mission integration. Not included here are gain reconstruction uncertainties, estimated to be 
of order 0.1 % . 

Table 6. Main characteristics of HFI full mission maps. 


Characteristic 


Reference frequency 

y [GHz] 



Notes 

100 

143 

217 

353 

545 

857 

Number of bolometers. 

8 

11 

12 

12 

3 

4 

al 

Effective beam FWHMi [arcmin] . 

9.68 

7.30 

5.02 

4.94 

4.83 

4.64 

bl 

Effective beam FWHM 2 [arcmin] . 

9.66 

7.22 

4.90 

4.92 

4.67 

4.22 

b2 

Effective beam ellipticity e . 

1.186 

1.040 

1.169 

1.166 

1.137 

1.336 

b3 

Noise per beam solid angle [//Kcmb]. 

7.5 

4.3 

8.7 

29.7 



cl 

[kJy sr"^]. 





9.1 

8.8 

cl 

Temperature noise fuKcMB deg]. 

1.29 

0.55 

0.78 

2.56 



c2 

[kJy sr"^ deg]. 





0.78 

0.72 

c2 

Polarization noise [//Kcmb deg]. 

1.96 

1.17 

1.75 

7.31 



c3 

Calibration accuracy [%] . 

0.09 

0.07 

0.16 

0.78 

l.l(+5) 

1.4(+5) 

d 

CIB monopole prediction [ MJy sr“^]. 

0.0030 

0.0079 

0.033 

0.13 

0.35 

0.64 

e 

Zodiacal light level correction [Kcmb]. 

4.3 X 10-’ 

9.4 X 10-^ 

3.8 X 10-6 3.4 X 10-5 



e2 

[MJysr'] _ 





0.04 

0.12 

e2 


Number of bolometers whose data were used in producing the channel map. 

FWHM of the Gaussian whose solid angle is equivalent to that of the effective beams. 

FWHM of the elliptical Gaussian fit. 

Ratio of the major to minor axis of the best-fit Gaussian averaged over the full sky. 

Estimate of the noise per beam solid angle, as given in bl. 

Estimate of the noise in intensity scaled to 1° assuming that the noise is white. 

Estimate of the noise in polarization scaled to 1° assuming that the noise is white. 

^ Calibration accuracy (at 545 and 857 GHz, the 5 % accounts for the model uncertainty). 

^ According to the Bethermin et al. (2012) model, whose uncertainty is estimated to be at the 20 % level (also for constant v/y). 
Zero-level correction to be applied on zodiacal-light corrected maps. 


and estimation of lensing potential (Planck Collaboration XV 
2015). We note here that the low-7’ polarization likelihood used 
in Planck Collaboration XIII (2015) is based exclusively on the 
70 GHz polarization map, cleaned of foregrounds by use of the 
30 and 353 GHz polarization maps. Both of these frequencies 
have been shown to be free of the kind of systematic errors 
that still affect intermediate frequencies (Planck Collaboration II 
2015; Planck Collaboration VIII 2015). 


S.2. CMB power spectra 

The foreground-subtracted, frequency-averaged, cross-half¬ 
mission TT spectrum is plotted in Fig. 11, together with the 
Commander power spectrum at multipoles i < 29. The figure 
also shows the best-fit base ACDM theoretical spectrum fitted to 
the PlanckTT-i-lowP likelihood, together with residuals (bottom 
panel) and +1 cr uncertainties. 
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Fig. 9. Maximum posterior CMB intensity map at 5' resolution derived from the joint baseline analysis of Planck, WMAP, and 
408 MHz observations. A small strip of the Galactic plane, 1.6 % of the sky, is filled in by a constrained realization that has the same 
statistical properties as the rest of the sky. 



-3 0 3 -3 0 3 

/iK /iK 


Fig. 10. Maximum posterior amplitude Stokes Q (left) and U {right) maps derived from Planck observations between 30 and 
353 GHz. These mapS have been highpass-filtered with a cosine-apodized filter between ^ = 20 and 40, and the a 17 % re¬ 
gion of the Galactic plane has been replaced with a constrained Gaussian realization (Planck Collaboration IX 2015). From 
Planck Collaboration X (2015). 


8.2.1. Polarization power spectra 

In addition to the TT spectra, the 2015 Planck likelihood in¬ 
cludes the TE and EE spectra. Figure 12 shows the TE and EE 
power spectra calculated from the 2015 data and including all 
frequency combinations. The theory curve shown in the figure 
is the best-fit base ACDM model fitted to the temperature spec¬ 


tra using the PlanckTT-l-lowP likelihood. The residuals shown 
in Fig. 12 are higher than expected and provide evidence of 
residual instrumental systematics in the TE and EE spectra. It 
is currently believed that the dominant source of errors is beam 
mismatch generating leakage from temperature to polarization 
at low levels of a few in We urge caution in the in- 
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terpretation of any features in these spectra, which should be 
viewed as work in progress. Nonetheless, we find a high level of 
consistency in results between the TT and the full TT+TE+EE 
likelihoods. Eurthermore, the cosmological parameters (which 
do not depend strongly on r) derived from the TE spectra have 
comparable errors to the TT-derived parameters, and they are 
consistent to within typically 0.5 cr or better. 

8.2.2. Number of modes 

One way of assessing the constraining power contained in a par¬ 
ticular measurement of CMB anisotropies is to determine the 
effective number of modes that have been measured. This 
is equivalent to estimating 2 times the square of the total S/N 
in the power spectra, a measure that contains all the available 
cosmological information if we assume that the anisotropies are 
purely Gaussian (and hence ignore all non-Gaussian informa¬ 
tion coming from lensing, the GIB, cross-correlations with other 
probes, etc.). Carrying out this procedure for the Planck 2013 
TT power spectrum data provided in Planck Collaboration XV 
(2014) and Planck Collaboration XVI (2014), yields the number 
826 000 (which includes the effects of instrumental noise, cos¬ 
mic variance and masking). The 2015 TT data have increased 
this value to 1 114 000, with TE and EE adding a further 60 000 
and 96 000 modes, respectively.^ Erom this perspective the 2015 
Planck data constrain approximately 55 % more modes than in 
the 2013 release. Of course this is not the whole story, since 
some pieces of information are more valuable than others, and 
in fact Planck is able to place considerably tighter constraints on 
particular parameters (e.g., reionization optical depth or certain 
extensions to the base ACDM model) by including new polar¬ 
ization data. 

8.2.3. Peaks in the power spectra 

The fidelity with which Planck has measured the and 

power spectra enables us to precisely estimate the underly¬ 
ing cosmological parameters (see Sect. 10), but the Qs are them¬ 
selves a set of cosmological observables, whose properties can 
be described independently of any model. The acoustic peaks in 
the Qs reveal the underlying physics of oscillating sound waves 
in the coupled photon-baryon fluid, driven by dark matter poten¬ 
tial perturbations, and one can talk about the fundamental mode, 
the first harmonic, and so on. Hence it is natural to ask about 
the positions of the individual peaks in the power spectra as em¬ 
pirical information that becomes part of the canon of facts now 
known about our Universe. 

Here we use the Planck data directly to fit for the multipoles 
of individual features in the measured TT, TE, and EE power 
spectra. We specifically use the CMB-only bandpowers given 
in Planck Collaboration XI (2015), adopting the same weight¬ 
ing scheme within each bin. Eitting for the positions and ampli¬ 
tudes of features in the bandpowers is a topic with a long history, 
with approaches becoming more sophisticated as the fidelity of 
the data improved (e.g., Scott & White 1994, Hancock & Rocha 
1997, Knox & Page 2000, de Bemardis et al. 2002, Bond et al. 
2003, Page et al. 2003, Durrer et al. 2003, Readhead et al. 2004, 
Jones et al. 2006, Hinshaw et al. 2007, Corasaniti & Melchiorri 
2008, Pryke et al. 2009). Eollowing earlier approaches, we fit 
Gaussians to the peaks in and but parabolas to the 

^Here we have used the basic (and conservative) likelihood; more 
modes are effectively probed by Planck if one includes larger sky frac¬ 
tions. 


peaks. We have to remove a featureless damping tail to 
fit the higher region and care has to be taken to treat 
the lowest-^ “recombination” peak in We explicitly fo¬ 
cus on peaks (ignoring the troughs) in the conventional quantity 
= l{l + l)Q/27r; note that other quantities (e.g., Q) will 
have maxima at slightly different multipoles, and that the choice 
of bandpowers to use for fitting each peak is somewhat subjec¬ 
tive. Our numerical values, presented in in Table 7, are consistent 
with previous estimates, but with a dramatically increased num¬ 
ber of peaks measured. Planck detects 19 peaks (with the eighth 
CY peak detection being marginal), and an essentially equiva¬ 
lent number of troughs. 

Table 7. Planck peak positions and amplitudes. 


Peak 

Number 

Position [€] 

Amplitude [yuK^] 

TT power spectrum 



First . 

220.0+ 0.5 

5717 + 35 

Second. 

537.5 + 0.7 

2582 + 11 

Third. 

810.8+ 0.7 

2523 + 10 

Fourth. 

1120.9+ 1.0 

1237 + 4 

Fifth. 

1444.2+ 1.1 

797.1 + 3.1 

Sixth. 

1776 + 5 

377.4+ 2.9 

Seventh . 

2081 +25 

214 + 4 

Eighth. 

2395 + 24 

105 + 4 

TE power spectrum 



First . 

308.5 + 0.4 

115.9 + 1.1 

Second. 

595.3 + 0.7 

28.6 + 1.1 

Third. 

916.9 + 0.5 

58.4+1.0 

Fourth. 

1224 + 1.0 

0.7 + 0.5 

Fifth. 

1536 +2.8 

5.6+ 1.3 

Sixth. 

1861 +4 

1.2 + 1.0 

EE power spectrum 



First . 

137 +6 

1.15 + 0.07 

Second. 

397.2 + 0.5 

22.04 + 0.14 

Third. 

690.8 + 0.6 

37.35 + 0.25 

Fourth. 

992.1 + 1.3 

41.8 +0.5 

Fifth. 

1296 + 4 

31.6 +1.0 


8.3. CMB lensing products 

Planck is the first experiment with the sky coverage, resolu¬ 
tion and sensitivity to form a full-sky reconstruction of the pro¬ 
jected mass, along every line of sight back to the surface of 
last scattering. Eigure 13 shows the 2015 Planck lensing map 
(Planck Collaboration XV 2015) using as input the CMB maps 
produced by the SMICA code. The map combines five possi¬ 
ble quadratic estimators based on the various correlations of the 
CMB temperature (T) and polarization {E and B). 

8.4. Likelihood code 
8.4.1. CMB likelihood 

We adopt the same general methodology for the 2015 likeli¬ 
hood as in 2013, extended to include Planck polarization data. 
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Fig. 11. The Planck 2015 temperature power spectrum. At multipoles ^ > 30 we show the maximum likelihood frequency-averaged 
temperature spectrum computed from the cross-half-mission likelihood with foreground and other nuisance parameters determined 
from the MCMC analysis of the base ACDM cosmology. In the multipole range 2 < ^ < 29, we plot the power spectrum estimates 
from the Coinmander component-separation algorithm computed over 94 % of the sky. The best-fit base ACDM theoretical spectrum 
fitted to the PlanckTT-i-lowP likelihood is plotted in the upper panel. Residuals with respect to this model are shown in the lower 
panel. The error bars show +1 cr uncertainties. From Planck Collaboration XIII (2015). 
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Fig. 12. Frequency-averaged TE (left) and EE (right) spectra (without fitting for T-P leakage). The theoretical TE and EE spectra 
plotted in the upper panel of each plot are computed from the best-fit model of Fig. 11. Residuals with respect to this theoretical 
model are shown in the lower panel in each plot. The error bars show +1 cr errors. The green lines in the lower panels show the 
best-fit temperature-to-polarization leakage model, fitted separately to the TE and EE spectra. From Planck Collaboration XIII 
(2015). 


The likelihood is a hybrid combination of a low-multipole pix- 
elbased likelihood with a high-multipole likelihood constructed 
from cross-spectra - see Planck Collaboration XI (2015) for de¬ 
tails. Note that we use the notation “Planck TT” when we are 
referring to the likelihood deriving from the TT spectrum. 

At low multipoles we now use Planck instead of WMAP for 
polarization information. The 70 GHz LFI polarization maps are 
cleaned with the LFI 30 GHz and HFI 353 GHz maps to mitigate 
foreground contamination. Based on null tests, the LFI-cleaned 


polarization maps are used over 46 % of the sky in the low mul¬ 
tipole likelihood (referred to as “lowP”). The Coinmander tem¬ 
perature solution, constructed from all Planck frequency maps, 
together with the Haslam 408 MHz and WMAP maps, is used 
over 93 % of the sky. The temperature and polarization data are 
then treated in a unified low-resolution pixel-based manner for 
the multipole range 7’ = 2 to 29. 

The high-^ likelihood uses pseudo-C^ cross-spectra from 
HFI 100, 143, and 217 GHz maps in a “fiducial Gaussian” ap- 
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Fig. 13. Wiener-filtered lensing potential estimate with minimal masking (using the NILC component separated map), in Galactic 
coordinates with a Mollweide projection (Planck Collaboration XV 2015). The reconstruction has been bandlimited to 8 < L < 2048 
(where, following convention, L is used as the multipole index in the lensing power spectrum). 


proximation, employing analytic covariance matrices calculated 
for a fiducial cosmological model. Unresolved foregrounds are 
modelled parametrically using power spectrum templates, with 
only minor changes to the model adopted in the 2013 analy¬ 
sis. The baseline high-multipole likelihood uses cross-spectra 
between frequency maps constructed from the first and second 
halves of the full mission data, to reduce any possible biases 
from co-temporal systematics. We also make more aggressive 
use of sky at all frequencies in the 2015 analysis. The most sig¬ 
nificant change is the addition of the option to include the TE 
and EE power spectra and the associated covariance matrices 
into the scheme, to form a combined TT, TE, EE likelihood at 
high multipoles (referred to as PlanckTT,TE,EE). Although we 
find firm evidence for systematics associated with temperature- 
to-polarization leakage in the TE and EE spectra, these system¬ 
atics are at low levels. We find a high level of consistency be¬ 
tween the TT, TE, and EE spectra for the cosmological models 
analyzed in the 2015 Planck papers. However, in this data re¬ 
lease, we regard the combined TT, TE, and EE Planck results 
as preliminary and hence recommend the TT likelihood as the 
baseline. 


8.4.2. Lensing likelihood 

Our power spectrum measurement constrains the lensing poten¬ 
tial power spectrum to a precision of +2.5 %, corresponding to a 
1.2 % constraint on the overall amplitude of matter fiuctuations 
(cTg), a measurement with considerable power for constraining 
cosmology. We have constructed two Gaussian bandpower like¬ 
lihoods based on the lensing power spectrum measurement de¬ 
scribed in Sect. 8.4.1 and plotted in Eig. 20. The first likeli¬ 
hood uses a conservative bandpower range, 40 < L < 400, 
with linear binning, following the temperature-only likelihood 
released in 2013. The second likelihood uses a more aggres¬ 
sive range with 8 < L < 2048, and bins that are linear in 


L^-^. Both likelihoods incorporate temperature and polarization 
data. We incorporate uncertainties in the estimator normaliza¬ 
tion and bias corrections directly into the likelihood, using pre¬ 
calculated derivatives of these terms with respect to the CMB 
temperature and polarization power spectra. The construction of 
the lensing likelihood is described in Planck Collaboration XV 
(2015), and its cosmological implications are discussed in detail 
in Planck Collaboration XIII (2015). 

9. Astrophysics products 

9.1. The Second Planck Catalogue of Compact Sources 

The Second Planck Catalogue of Compact Sources (PCCS2; 
Planck Collaboration XXVI 2015) is the catalogue of sources 
detected from the full duration of Planck operations, referred to 
as the “extended” mission. It consists of compact sources, both 
Galactic and extragalactic, detected over the entire sky. Compact 
sources are detected in the single-frequency maps and assigned 
to one of two sub-catalogues, the PCCS2 or PCCS2E. The first 
of these allows the user to produce additional sub-catalogues at 
higher reliabilities than the target 80 % reliability of the full cat¬ 
alogue. The second list contains sources whose reliability cannot 
be estimated, because they are embedded in a bright and com¬ 
plex (e.g. filamentary) background of emission. 

The total number of sources in the catalogue ranges from 
1560 at 30 GHz up to 48181 sources at 857 GHz. Both sub¬ 
catalogues include polarization measurements, in the form of 
polarized fiux densities and orientation angles, or upper-limits, 
for all seven polarization-sensitive Planck channels. The num¬ 
ber of sources with polarization information (other than upper- 
limits) in the catalogue ranges from 113 in the lowest polar¬ 
ized frequency channel (30 GHz) up to 666 in the highest po¬ 
larized frequency channel (353 GHz). The improved data pro¬ 
cessing of the full-mission maps and their reduced instrumen¬ 
tal noise levels allow us to increase the number of objects in 
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the catalogue, improving its completeness for the target 80 % 
reliability as compared with the previous versions, the PCCS 
(Planck Collaboration XXVIII 2014) and the Early Release 
Compact Source Catalogue (ERCSC; Planck Collaboration XIII 
2011). The improvements are most pronounced for the LEI chan¬ 
nels, due to the much larger increase in the data available. The 
completeness of the 857 GHz channel, however, has not im¬ 
proved; this is due to a more refined reliability assessment, which 
resulted in a higher S/N threshold being applied in the selection 
function of this catalogue. The reliability of the PCCS2 cata¬ 
logue at 857 GHz, however, is higher than that of the PCCS. 

9.2. The Second Planck Catalogue of Clusters 

The Second Planck Catalogue of SZ Sources (PSZ2; 
Planck Collaboration XXVII 2015), based on the full mission 
data, is the largest SZ-selected sample of galaxy clusters yet pro¬ 
duced and the deepest all-sky catalogue of galaxy clusters. It 
contains 1653 detections, of which 1203 are confirmed clusters 
with identified counterparts in external data sets, and is the first 
SZ-selected cluster survey containing > 10^ confirmed clusters. 
A total of 937 sources from the half-mission catalogue (PSZl) 
released in 2013 are included, as well as 716 new detections. 
The completeness, which is provided as a product with the cat¬ 
alogue, is determined using simulated signal injection, validated 
through comparison to external data, and is shown to be consis¬ 
tent with semi-analytic expectations. The reliability is character¬ 
ized using high-fidelity simulated observations and a machine¬ 
learning-based quality assessment, which together place a robust 
lower limit of 83 % on the purity. Using simulations, we find 
that the F 5 R 500 estimates are robust to pressure-profile variations 
and beam systematics; however, accurate conversion to 1500 re¬ 
quires the use of prior information on the cluster extent. Results 
of a multi-wavelength search for counterparts in ancillary data 
(which makes use of radio, microwave, infra-red, optical, and 
X-ray data-sets, and which places emphasis on the robustness of 
the counterpart match) are included in the catalogue. We discuss 
the physical properties of the new sample and identify a popula¬ 
tion of low-redshift X-ray under-luminous clusters revealed by 
SZ selection. Eigure 14 shows the masses and redshifts for the 
1093 PSZ2 clusters with known redshifts. 

9.3. The Planck Catalogue of Galactic Cold Clumps 

The Planck catalogue of Galactic Cold Clumps (PGCC, 
Planck Collaboration XXVIII 2015) contains Galactic sources 
that have been identified as cold sources in Planck data. We 
ran the CoCoCoDeT (Montier et al. 2010) multi-frequency point 
source detection algorithm on the Planck 857, 545, and 353 GHz 
data and the IRIS 3 THz data (Miville-Deschenes & Lagache 
2005), at a resolution of 5'. This selects point sources exhibit¬ 
ing submillimetre excess in the 857, 545, and 353 GHz Planck 
bands simultaneously, compared to the average colour of the 
background, which is typical of sources appearing colder than 
their environment. 

The PGCC catalogue is the full version of the Early Cold 
Core (ECC) catalogue released in 2011, which was part of 
the ERCSC (Planck Collaboration VII 2011). The ECC cata¬ 
logue was built on the first 295 days of Planck data, and con¬ 
tains 915 sources selected to ensure T < 14K and S/N> 15. 
A statistical description of the ECC and the extended cata¬ 
logue (including sources at all temperatures and with S/N > 4) 
is given in Planck Collaboration XXII (2011), while a detailed 



Redshift 


Fig. 14. Distribution of the 1093 PSZ2 clusters with counterparts 
with known redshift in the M 500 -Z plane. New PSZ2-detected 
clusters are indicated with red dots, while commmon PSZl and 
PSZ2 clusters are indicated by black dots. Green dots mark the 
common PSZ2-PSZ1 detections where we have updated the red¬ 
shift value for the PSZ2. The solid, dashed, and dotted lines in¬ 
dicate the limiting mass at 20 %, 50 %, and 80 % survey com¬ 
pleteness, respectively. 


description of a subsample of 10 sources was presented in 
Planck Collaboration XXIII (2011). The PGCC catalogue, in¬ 
cluded in the 2015 Planck release, has now been built on the full 
Planck mission data, and contains 13 188 Galactic sources, plus 
54 sources located in the Large and Small Magellanic Clouds. 

The morphology of each source is obtained using a Gaussian 
elliptical fit, which is then used to estimate fiux densities in 
all bands through aperture photometry. Depending on the S/N 
of the fiux density estimates, three categories of sources are 
identified: 6993 sources with reliable fiux densities in all bands 
(FLUX_QUALITY=1); 3755 sources with fiux density estimates in 
all bands except 3 THz (FLUX_QUALITY=2), which are consid¬ 
ered very cold candidates; and 2440 sources without reliable fiux 
density estimates (FLUX_QUALITY=3), usually due to a complex 
environment, which are considered poor candidates. 

Distance estimates have been obtained for 5574 PGCC 
sources by combining seven different methods. While PGCC 
sources are mainly located in the solar neighbourhood, with 
88 % of sources with reliable distance estimates lying within 
2 kpc of the Sun, distance estimates range from a few hundred 
parsecs towards local molecular clouds to 10.5 kpc towards the 
Galactic centre. 

The temperature of each source is obtained by fitting a mod¬ 
ified blackbody to the spectral energy density from 3 THz to 
353 GHz, considering the spectral index yS as a free parameter 
when possible. PGCC sources have an average temperature of 
13-14.5 K, depending on fiux quality category, and range from 
5.8 to 20 K. Other physical parameters have been derived, such 
as the H 2 column density, the physical size, the mass, the den¬ 
sity, and the luminosity. It appears that the PGCC contains a 
large variety of objects with very different properties, from com¬ 
pact and dense cores to large and massive molecular clouds, lo¬ 
cated all over the sky. While a large Herschel program (HKP- 
GCC) already followed-up 315 PGCC sources with the PACS 
and SPIRE instruments, the PGCC catalogue is the first all-sky 
sample of Galactic cold sources obtained with a homogeneous 
method, and hence represents a goldmine for investigations of 
the early phases of star formation in various environments. 
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Fig. 15. All-sky distribution of the 13 188 PGCC Galactic cold 
clumps (black dots) and the 54 cold sources (grey dots) located 
in the Large and Small Magellanic Clouds. The background map 
is the 857 GHz Planck band, shown in logarithmic scale from 
10-2 to lO^MJysr-i. 


9.4. Diffuse Galactic foregrounds from CMB component 
separation 

As was done in 2013, in Planck Collaboration X (2015) we es¬ 
tablish a single parametric model of the microwave sky, ac¬ 
counting simultaneously for all significant diffuse astrophys- 
ical components and relevant instrumental effects using the 
Bayesian Coimnander analysis framework (Eriksen et al. 2004, 
2006, 2008). The 2015 analysis is extended in multiple direc¬ 
tions. First, instead of 15.5 months of temperature data, the new 
analysis includes the full Planck mission data—50 months of 
LFI and 29 months of HFI data—in both temperature and po¬ 
larization. Second, we now also include the 9-year WMAP ob¬ 
servations between 23 and 94 GHz (Bennett et al. 2013) and a 
408 MHz survey map (Haslam et al. 1982), providing enough 
frequency constraints to decompose the low-frequency fore¬ 
grounds into separate synchrotron, free-free, and spinning dust 
components. Third, we now include the Planck 545 and 857 GHz 
frequency bands, allowing us to constrain the thermal dust tem¬ 
perature and emissivity index with greater precision, thereby 
reducing degeneracies between CMB, CO, and free-free emis¬ 
sion. Fourth, the present analysis implements a multi-resolution 
strategy to provide component maps at high angular resolu¬ 
tion. Specifically, the CMB is recovered with angular resolution 
5' FWHM (Planck Collaboration IX 2015), thermal dust emis¬ 
sion and CO 7 = 2^1 lines are recovered at 7.'5 FWHM, and 
synchrotron, free-free, and spinning dust are recovered at 1° 
FWHM. 

An important difference with respect to 2013 is that we em¬ 
ploy individual detector and detector set maps as inputs, instead 
of fully combined frequency maps. The increase in the number 
of input maps allows us to make many null tests that are used to 
reject individual maps exhibiting significant levels of systematic 
effects. In addition, in our analysis we allow our model to fit for 
two important instrumental effects: relative calibration between 
detectors; and bandpass uncertainties. 

The sum of these improvements allows us to reconstruct a 
total of six primary emission mechanisms in temperature: CMB; 
synchrotron; free-free; spinning dust; CO; and thermal dust 
emission—in addition to two secondary components, namely 
thermal SZ emission around the Coma and Virgo regions, and 
molecular line emission between 90 and 100 GHz. For polariza¬ 
tion, we reconstruct three primary emission mechanisms: CMB; 
synchrotron; and thermal dust. All of these components are de¬ 
livered as part of the 2015 Planck release. 


Figures 16 and 17 (Planck Collaboration X 2015) show the 
diffuse high-latitude Galactic foreground components deter¬ 
mined from component separation in temperature and polar¬ 
ization. Figure 18 shows the spectra of fluctuations of diffuse 
foreground components in temperature and polarization, com¬ 
pared to the CMB. The sky model presented in this paper pro¬ 
vides an impressive fit to the current data, with temperature 
residuals at the few microkelvin level at high latitudes across 
the CMB-dominated frequencies, and with median fractional er¬ 
rors below 1 % in the Galactic plane across the Planck frequen¬ 
cies. For polarization, the residuals are statistically consistent 
with instrumental noise at high latitudes, but limited by signif¬ 
icant temperature-to-polarization leakage in the Galactic plane. 
Overall, this model represents the most accurate and complete 
description currently available of the astrophysical sky between 
20 and 857 GHz. 


9.5. Carbon monoxide emission 

Carbon monoxide emission lines are present in all HFI fre¬ 
quency bands, except 143 GHz. Using component separa¬ 
tion techniques, the three lowest rotational transitions can 
be extracted from Planck data, providing full sky maps 
of the CO7=1^0, 7 = 2^1, and 7 = 3^2 transitions 
(Planck Collaboration XIII 2014). For the 2015 release, data 
from the full mission and better control of systematic errors 
lead to better maps. Table 8 summarizes the products. Figure 16 
shows the Commander maps of all three transitions. 


- Type 1 maps are produced by a single-channel analysis, 
where individual bolometer maps are linearly combined to 
produce maps of the CO(l^O), CO(2^1), and CO(3^2) 
emission lines at the native resolution of the Planck maps. 
Although noisier than the other approaches, using informa¬ 
tion from a single channel strongly limits contamination 
from other Galactic components, such as dust or free-free 
emission. This makes Type 1 maps suitable for studying 
emission in the Galactic disk and CO-rich regions, but not 
for the high-Galactic latitudes where the CO emission is be¬ 
low the noise level. 

- Type 2 maps of CO(l^O) and CO(2^1) have been obtained 
using multi-channel information (i.e., using linear combina¬ 
tion of Planck channel maps smoothed to 15'). Using fre¬ 
quency maps, this type of map has a higher signal-to-noise 
ratio, allowing for their use in fainter high-Galactic latitude 
regions. They are, however, more susceptible to dust con¬ 
tamination, especially for CO(2^1), which makes them less 
suitable in the Galactic plane than Type 1 maps. 

- A high-resolution Type 3 map, as defined in 
Planck Collaboration XIII (2014), is not being deliv¬ 
ered in the 2015 data release. Alternatively, another set of 
CO maps has been produced as part of the full Coimnander 
baseline multi-component model, which is described in 
Planck Collaboration X (2015). 


Type 1 and Type 2 maps are released with associated stan¬ 
dard deviation maps, error maps, and masks. The suite of tests 
detailed in Planck Collaboration XIII (2014) has been repeated 
on the new Type 1 and Type 2 maps, which have been found to 
perform as well as their 2013 counterparts, even though small 
variations (^ 2 to 5 Kkm s“^) exist in the Galactic plane. 
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Fig. 16. Maximum posterior intensity maps derived from the joint analysis of Planck, WMAP, and 408 MHz observations 
(Planck Collaboration X 2015). From left to right, top to bottom: CMB; synchrotron; free-free; spinning dust; thermal dust; line 
emission around 90 GHz; CO / = 1 ^ 0; CO / = 2 ^ 1, and CO J = 3^2. 

Table 8. Summary of main CO product characteristics. 


Map 

Algorithm 

CO line 

Resolution 

[arcmin] 

Noise rms [Krj km s 

Analysis details 

15' FWHM 60' FWHM 

Frequencies [GHz] 

Model 

Type 1 . . . . 

MILCA 

/=1^0 

9.6 

1.4 

0.34 

100 (bol maps) 

CO, CMB 


MILCA 

7 = 2^1 

5.0 

0.53 

0.16 

217 (bol maps) 

CO, CMB, dust 


MILCA 

7 = 3^2 

4.8 

0.55 

0.18 

353 (bol maps) 

CO, dust 

Type 2 . . . . 

MILCA 

7=1^0 

15 

0.39 

0.085 

70, 100, 143, 353 

CO, CMB, dust, free-free 


MILCA 

7 = 2^1 

15 

0.11 

0.042 

70, 143, 217, 353 

CO, CMB, dust, free-free 


Commander 

7=1^0 

60 


0.084 

0.408-857 

Full 


Commander 

7 = 2^1 

60 


0.037 

0.408-857 

Full 


Commander 

7 = 3^2 

60 


0.060 

0.408-857 

Full 

Type 3 ... . 

Commander 

7 = 2^F 

7.5 

0.090 

0.031 

143-857 

CO, CMB, dust 


Commander-Ruler 

7=1^0*^’" 

5.5 

0.19 

0.082 

30-353 

CO, CMB, dust, low-freq 


^ Formally a weighted average of CO / = 2 ^ 1 and J -3^2, but strongly dominated by CO / = 2 ^ 1. 

^ Formally a weighted average of CO / =1^0, / = 2^1 and J -3^2, but strongly dominated by CO / = 1 ^ 0. 
^ Only published in 2013. 


9.5.1. All-sky Sunyaev-Zeldovich emission 

The 30 to 857 GHz frequency channel maps from the Planck 
satellite survey were used to construct an all-sky map of the ther¬ 
mal Sunyaev-Zeldovich (tSZ) effect Planck Collaboration XXII 
(2015). As discussed in Planck Collaboration XXI (2014), 
we apply to those maps specifically tailored component 
separation algorithms, MILCA Hurier et al. (2010) and NILC 


Remazeilles et al. (2011), that allow us to separate the tSZ 
emission from foreground contamination, including the CMB. 
An orthographic view of the reconstructed Compton y-map 
in Healpix pixelization is presented in Fig. 19. This y- 
map has been characterized in terms of noise properties and 
residual foreground contamination, mainly thermal dust emis¬ 
sion at large angular scales and CIB and extragalactic point 
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Fig. 17. Maximum posterior amplitude polarization maps derived from the Planck observations between 30 and 353 GHz 
(Planck Collaboration X 2015). The left and right columns show the Stokes Q and U parameters, respectively. Rows show, from top 
to bottom: CMB; synchrotron polarization at 30 GHz; and thermal dust polarization at 353 GHz. The CMB map has been highpass- 
filtered with a cosine-apodized filter between ^ = 20 and 40, and the Galactic plane (defined by the 17 % CPM83 mask) has been 
replaced with a constrained Gaussian realization (Planck Collaboration IX 2015). 



30 44 70 100 143 217 353 



Frequency (GHz) 


Fig. 18. Brightness temperature rms of the high-latitude sky as a function of frequency and astrophysical component for temperature 
(left) and polarization {right). For temperature, each component is smoothed to an angular resolution of 1° FWHM, and the lower 
and upper edges of each line are defined by masks covering 81 and 93 % of the sky, respectively. For polarization, the corresponding 
smoothing scale is 40', and the sky fractions are 73 and 93 %. 


sources at small angular scales. Blindly-detected clusters in 
this map are consistent with those from the PSZ2 catalogue 
Planck Collaboration XXVII (2015), both in terms of cluster 
numbers and integrated flux. Furthermore, by stacking individu¬ 
ally undetected groups and clusters of galaxies we find that the 
y-map is consistent with tSZ emission even for low S/N regions. 
Using foreground models derived in Planck Collaboration XXIII 
(2015) we are able to measure the tSZ angular power spectrum 
over 50 % of the sky. We conclude that the y-map is dominated 
by tSZ signal in the multipole range, 20 < £ < 800. Similar re¬ 


sults are obtained from a high-order statistic analysis. The recon¬ 
structed y-map is delivered as part of the Planck 2015 release. 
We also deliver a foreground mask (with known point sources 
and regions with strong contamination from Galactic emission 
masked out), a noise variance map, the estimated power spec¬ 
trum, and the weights for the NILC algorithm. 
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Fig. 19. Orthographic projection of the reconstructed Planck all-sky y-map in Compton parameter units (Planck Collaboration XXII 
2015). For illustration purposes and to enhance the tSZ signal-to-noise ratio, the y-map has been Wiener filtered. Positive sources in 
the map correspond to clusters and super-clusters of galaxies with strong tSZ emission. In particular, the Coma and Virgo clusters 
are clearly visible near the north Galactic pole. The region of strongest contamination from Galactic foreground emission in the 
Galactic plane has been partially masked. 


10. Planck 2015 cosmology results 

Since their discovery, anisotropies in the CMB have contributed 
significantly to defining our cosmological model and measuring 
its key parameters. The standard model of cosmology is based 
upon a spatially fiat, expanding Universe whose dynamics are 
governed by General Relativity and dominated by cold dark mat¬ 
ter and a cosmological constant (A). The seeds of structure have 
Gaussian statistics and form an almost scale-invariant spectrum 
of adiabatic fluctuations. The 2015 Planck data remain in excel¬ 
lent agreement with this paradigm, and continue to tighten the 
constraints on deviations and reduce the uncertainty on the key 
cosmological parameters. 

The major methodological changes in the steps going 
from sky maps to cosmological parameters are discussed 
in Planck Collaboration XII (2015); Planck Collaboration XIII 
(2015). These include the use of Planck polarization data in¬ 
stead of WMAP, changes to the foreground masks to include 
more sky and dramatically reduce the number of point source 
“holes,” minor changes to the foreground models, improve¬ 
ments to the data processing and use of cross-half-mission likeli¬ 
hoods (Planck Collaboration XI2015; Planck Collaboration XII 
2015). We find good agreement with our earlier results, with in¬ 
creased precision. 

10.1. Cosmological parameters 

Planck's measurements of the cosmological parameters de¬ 
rived from the full mission are presented and discussed in 
Planck Collaboration XIII (2015). As in our previous release, 
the data are in excellent agreement with the predictions of 


Table 9. Parameter 68 % confidence levels for the base ACDM 
cosmology computed from the Planck CMB power spectra, in 
combination with the CMB lensing likelihood (“lensing”). 


Parameter P/a^zck TT+lowP+lensing 


. 0.02226 + 0.00023 

. 0.1186 + 0.0020 

lOO^Mc. 1.04103 + 0.00046 

T. 0.066 + 0.016 

ln(10^%). 3.062 + 0.029 

n, . 0.9677 + 0.0060 

Ho . 67.8 + 0.9 

. 0.308 + 0.012 

. 0.1415 + 0.0019 

. 0.09591 + 0.00045 

0-8. 0.815 + 0.009 

. 0.4521 + 0.0088 

Age/Gyr . 13.799 + 0.038 

r^rag. 147.60 + 0.43 

y^eq . 0.01027 + 0.00014 


the 6-parameter ACDM model (see Table 9), with parame¬ 
ters tightly constrained by the angular power spectrum. The 
best-fit model parameters from the full mission are typically 
within a fraction of a standard deviation of their results from 
Planck Collaboration XVI (2014), with no outliers. The con¬ 
straints on the parameters of the base ACDM model have im¬ 
proved by up to a factor of 3. The largest shifts are in the 
scalar spectral index, n^, which has increased by 0.7 cr, and the 
baryon density, tUb, which has increased by 0.6 cr. Both of these 


27 






















































Planck Collaboration: The Planck mission 


shifts are partly due to correction of a systematic error that con¬ 
tributed to a loss of power near I = 1800 in the 2013 results 
Planck Collaboration XIII (2015). This systematic also biased 
the inferences on Hq slightly low (by less than 0.5 cr). In addi¬ 
tion, the overall amplitude of the observed spectrum has shifted 
upwards by 2 % (in power), due to a calibration change, and 
the optical depth to Thomson scattering, r, has shifted down by 
nearly 1 cr. These shifts approximately cancel in the derived nor¬ 
malization of the matter power spectrum. The remaining shifts 
are consistent with the known changes in noise level, time- 
stream filtering, absolute calibration, beams, and other aspects 
of the data processing. 

Both the angular size of the sound horizon, 6^, and the cold 
dark matter density, tUc, have become significantly better deter¬ 
mined. The data at high-^ are now so precise, and the polariza¬ 
tion data so constraining, that we not only see very strong evi¬ 
dence for three species of light neutrinos, but can measure the ef¬ 
fective viscosity of the neutrino “fluid” to be non-zero at the 9 cr 
level. The constraint on the baryon density, tUb, is now compa¬ 
rable with the best quoted errors from big bang nucleosynthesis 
and suggests the possibility of calibrating nuclear capture cross- 
sections from CMB observations. The addition of polarization 
data has improved by an order of magnitude our upper limit on 
the annihilation rate of dark matter. 

Despite trying a wide range of extensions to the basic, 6 - 
parameter ACDM model, we find no significant evidence for 
a failure of the model. Within each extension of the parame¬ 
ter space, the default parameter values for the ACDM model 
remain a good fit to the data. This continues to hold when we 
combine the Planck data with other measurements, such as the 
distance scale measured by baryon acoustic oscillations (BAO) 
in galaxy surveys or Type la supernovae, or the growth of struc¬ 
ture determined by redshift-space distortions. Since our best- 
fit cosmology has shifted by very little since our 2013 release, 
we continue to see tensions with some analyses of other as- 
trophysical data sets (e.g., the abundance of clusters of galax¬ 
ies and weak gravitational lensing of galaxies or cosmic shear, 
and distances measured by BAO in the Lya forest at high-z). 
Planck Collaboration XIII (2015) shows that these tensions can¬ 
not be resolved with standard single parameter extensions of the 
base ACDM model. Resolving these discrepancies remains an 
area of active research. 

10.2. Constraints from large angular scales 

The anisotropy at large angular scales, particularly the polariza¬ 
tion, allows us to place tight constraints on the optical depth 
to Thomson scattering, r, and the tensor-to-scalar ratio, r. The 
Planck temperature data, in combination with CMB lensing and 
low-^ polarization measured at 70 GHz, prefer a lower opti¬ 
cal depth, T = 0.066 + 0.016, than the earlier inference from 
WMAP9 (r 0.09, which was used in our 2013 analysis), which 
implies a lower redshift of reionization (zre = 8 . 8 ^} 4 ). When 
cleaned of foregrounds using our 353 GHz channel, the WMAP 
polarization data are in good agreement with a lower optical 
depth. With the dramatic improvement in our CMB lensing de¬ 
tection, we are able to independently constrain r, finding compa¬ 
rably tight and consistent results (r = 0.071 + 0.016) without the 
use of low-^ polarization. This provides additional confidence in 
the results. 

While improved constraints on polarization at low-^ will 
eventually allow us to study the reionization epoch in more de¬ 
tail, at present the largest impact of the change in r comes from 
the implied downward shift in the inferred matter power spec¬ 


trum normalization, erg. As it happens, much of the downward 
shift in this parameter is largely cancelled by the upward shift in 
the CMB spectrum arising from the improved calibration in the 
current data release. 

Gravitational waves entering the horizon between recombi¬ 
nation and today give a “tensor” contribution to the large-scale 
temperature and polarization anisotropies. Our strongest Planck- 
only constraint still comes from temperature anisotropies at 
^ < 10^ (or k ^ 0.01 Mpc“^), and is thus limited by cos¬ 
mic variance and model-dependent. Tensor modes also gener¬ 
ate a 5-mode signal, which peaks at 7’ 10^, slightly smaller 

scales than the bulk of the temperature signal. The cosmolog¬ 
ical landscape became more complicated earlier this year with 
the detection of 5-mode polarization anisotropy by the BICEP2 
team (Ade et al. 2014). Analysis of Planck polarization data 
at high Galactic latitudes demonstrated that no region of the 
sky can be considered dust-free when searching for primor¬ 
dial 5-modes (Planck Collaboration Int. XXX 2014), and a joint 
analysis of BICEP2/Keck Array observations and Planck po¬ 
larization data (BICEP2/Keck Array and Planck Collaborations 
2015) shows that polarized dust emission contributes a signifi¬ 
cant part of the BICEP2 signal. Combining the Planck and re¬ 
vised BICEP2/Keck Array likelihoods leads to a 95 % upper 
limit of ro.oo 2 < 0.09. This eliminates any tension between the 
BICEP2 and Planck results, and in combination with our other 
constraints disfavours inflationary models with a 0 ^ potential. 
This and other implications for infiationary models in the early 
Universe are discussed more fully in Planck Collaboration XIII 
(2015); Planck Collaboration XX (2015). 

10.3. Dark energy and modified gravity 

Even though much of the weight in the Planck data lies at high 
redshift, Planck can still provide tight constraints on dark en¬ 
ergy and modified gravity, especially when used in combination 
with other probes. This is explored in Planck Collaboration XIV 
(2015), which focuses on tests of dark energy and modi¬ 
fied gravity on the scales where linear theory is most appli¬ 
cable, since these are the most theoretically robust. As for 
Planck Collaboration XIII (2015), the results are consistent with 
the simplest scenario, ACDM, though all constraints on dark 
energy models (including minimally-coupled scalar field mod¬ 
els or evolving equation of state models) and modified grav¬ 
ity models (including effective field theory, phenomenological, 
/(5), and coupled dark energy models) are considerably im¬ 
proved with respect to past analyses. In particular, we improve 
significantly the constraint on the density of dark energy at 
early times, finding that it has to be below 2 % (95 % confi¬ 
dence) of the critical density, even if it only plays a role below 
z = 50. Constraints are tighter if early dark energy is present 
since recombination, with Qg < 0.0071 (for the data combi¬ 
nation PlanckTT-rlensing-hBAO-rSNe-rTfo), and an even tighter 
bound if high-7 polarization is included. In models where per¬ 
turbations are modified, even if the background is ACDM, a few 
tensions appear, mainly driven by external data sets. 

10.4. Lensing of the CMB 

The CMB fluctuations measured by Planck provide a slightly 
perturbed image of the last-scattering surface, due to the effects 
of gravitational lensing by large-scale structure. Lensing slightly 
washes out the acoustic peaks of the CMB power spectrum, an 
effect we see in the Planck data at high significance. Lensing 
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also introduces distinctive non-Gaussian features into the CMB 
maps, which allow us to map and make statistical measure¬ 
ments of the gravitational potentials, and the associated matter. 
These are studied in detail in Planck Collaboration XV (2015). 
The lensing signal is consistent with the basic, 6 -parameter, 
ACDM model that best fits the temperature data. This gives us 
a very strong consistency check on the gravitational instabil¬ 
ity paradigm and the growth of structure over more than two 
decades in expansion factor. 

Since it provides sensitivity to the growth of structure be¬ 
tween the surface of last scattering and the present epoch, the 
lensing signal allows us to measure a number of important pa¬ 
rameters by breaking parameter degeneracies. Figure 20 shows 
the lensing power spectrum, which for the first time is measured 
with higher accuracy than it is predicted by the base ACDM 
model that fits the temperature data. With the temperature-only 
nominal mission data from the 2013 Planck data release, we 
were able to make the most powerful measurement of lensing 
to that date (at a level of 25 cr). In the current release, incorpo¬ 
rating additional temperature data, as well as entirely new po¬ 
larization information, we have nearly doubled the power of this 
measurement to 40 cr. This is the most significant detection to 
date, allowing lensing to be used as part of our precision cos¬ 
mology suite. 


10.5. Inflation 

The release of the 2013 Planck data and findings had an enor¬ 
mous impact upon the infiationary community, and the Planck 
2015 results continue to stress the importance of this win¬ 
dow into the early Universe. Planck Collaboration XX (2015) 
presents our constraints on inflationary models. The Planck data 
are consistent with a purely adiabatic, power-law spectrum of 
initial fluctuations, whose spectral index {n^ = 0.9677 + 0.006) 
is significantly different from unity. The addition of polarization 
data has significantly improved the constraints on any isocur¬ 
vature modes, which are now constrained at the percent level. 
Despite a detailed search, and study of several models, we see 
no statistically significant evidence for departures from a power 
law. The combination of Planck data with the BICEP2/Keck 
Array data provide a strong upper limit on the tensor-to-scalar 
ratio, and disfavour all monomial models (y(0) oc 0 ^^) with 
p>\. This is an important milestone, since these form the sim¬ 
plest class of inflationary models. 


10.6. Primordial non-Gaussianity 

(Planck Collaboration XVII 2015) for the first time uses polar¬ 
ization information to constrain non-Gaussian signals left by 
primordial physics. The results significantly reduce the allowed 
model space spanned by local, equilateral, and orthogonal non- 
Gaussianity, tighenting constraints by up to 45 %. In particular, 
= 0.8 ± 5.0, = -4 ± 43, and = -26 ± 21. 

In addition, the Planck 2015 analysis covers a greatly extended 
range of primordial 3-point and 4-point signals, constraining in¬ 
flationary model space as well as some proposed alternatives to 
inflation. The global picture that emerges is one of consistency 
with the premises of ACDM cosmology, namely that the struc¬ 
ture we observe today is the consequence of the passive evo¬ 
lution of adiabatic, Gaussian, nearly scale-invariant, primordial 
seed perturbations. 


10.7. Isotropy and statistics 

The Planck 2013 results determined the presence of statisti¬ 
cally anisotropic signals in the CMB, confirming previous stud¬ 
ies made using WMAP data. Such anomalies therefore consti¬ 
tute real features of the microwave sky, and potentially chal¬ 
lenge fundamental assumptions of the standard cosmological 
model. Planck Collaboration XVI (2015) extends these studies 
based mainly on the full Planck mission for temperature, but 
also including some polarization measurements. A large number 
of statistical tests indicate consistency with Gaussianity, while 
a power deficit at large angular scales is manifested in several 
ways, for example low map variance. The well-known ‘Cold 
Spot’ is identified through various methods. Tests of directional¬ 
ity suggest the presence of angular clustering from large to small 
scales, but at a significance that is dependent on the details of 
the approach. On large-angular scales, a dipolar power asymme¬ 
try is investigated through several approaches, and we address 
the subject of a posteriori correction. Our ability to include re¬ 
sults based on polarization data is limited by two factors. First, 
CMB polarization maps have been high-pass filtered to miti¬ 
gate residual large-scale systematic errors in the HFI channels, 
thus eliminating structure in the maps on angular scales larger 
than about 10°. Second, an observed noise mismatch between 
the simulations and the data prevents robust conclusions from 
being reached based on the nulll-hypothesis approach adopted 
throughout the paper. Nevertheless, we perform the first exami¬ 
nation of polarization data via a stacking analysis, in which the 
stacking of the data themselves necessarily acts to lower the ef¬ 
fect of the noise mismatch. We find that the morphology of the 
stacked peaks is consistent with the expectations of statistically 
isotropic simulations. Further studies of the large angular scale 
structure of the CMB polarization anisotropy will be conducted 
with data of improved quality expected to be released in 2016. 

10.8. Cosmology from clusters 

In 2013 we found an apparent tension between our primary CMB 
constraints and those from the Planck cluster counts, with the 
clusters preferring a lower normalization of the matter power 
spectrum, erg. The comparison is interesting because the cluster 
counts directly measure erg at low redshift and hence any tension 
could signal the need for extensions of the base model, such as 
non-minimal neutrino mass. However, limited knowledge of the 
normalization of the scaling relation between SZ signal and mass 
(usually called “mass bias”) continues to hamper the interpreta¬ 
tion of this result. 

Our 2015 cluster analysis benefits from a larger catalogue 
(438 objects versus the 189 in 2013), greater control of the se¬ 
lection function, and recent gravitational lensing determinations 
of the mass bias for Planck clusters. With the larger sample, 
we now fit the counts in the two-dimensional plane of redshift 
and S/N, allowing us to simultaneously constrain the slope of 
the scaling relation and the cosmological parameters. We exam¬ 
ine three new empirical determinations of the mass bias from 
gravitational lensing: Weighing the Giants (von der Linden et al. 
2014, WtC); the Canadian Cluster Comparison Project (CCCP; 
Hoekstra et al., private communication); and results from a new 
method based on CMB lensing (Melin & Bartlett 2014). We use 
these three results as priors because they measure the mass scale 
directly on samples of Planck clusters. 

The cluster constraints on erg and Qm are statistically iden¬ 
tical to those of 2013 when adopting the same scaling relation 
and mass bias; in this sense, we confirm the 2013 results with 
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Fig. 20. Lensing potential power spectrum estimate from the 2015 data release (Planck Collaboration XV 2015), based on the SMICA 
CMB map, as well as previous reconstructions from Planck as well as other experiments for comparison. 


the larger 2015 catalogue. Applying the three new mass bias pri¬ 
ors, we find that the WtG calibration reduces the tension with the 
primary CMB constraints to slightly more than 1 cr in the base 
model, and CCCP results in tension at just over 2 cr, similar to 
the case for the CMB lensing calibration. More detailed discus¬ 
sion of constraints from Planck cluster counts can be found in 
Planck Collaboration XXIV (2015). 


11. Planck 2015 astrophysics results 

11.1. Low frequency foregrounds 

Planck Collaboration XXV (2015) discusses Galactic 
foreground emission between 20 and 100 GHz, based 
primarily on the Commander component separation of 
(Planck Collaboration X 2015). The total intensity in this 
part of the spectrum is dominated by free-free and spinning 
dust emission, while polarization is dominated by synchrotron 
emission. 

Comparison with radio recombination line templates verifies 
the recovery of the free-free emission along the Galactic plane. 
Comparison of the high-latitude Ha emission with our free-free 
map shows residuals that correlate with dust optical depth, con¬ 
sistent with a fraction (~30 %) of Ha having been scattered by 
high-latitude dust. We highlight a number of diffuse spinning 
dust morphological features at high latitude. There is substantial 
spatial variation in the spinning dust spectrum, with the emission 
peak (in ly) ranging from below 20 GHz to more than 50 GHz. 
There is a strong tendency for the spinning dust component near 


many prominent H ii regions to have a higher peak frequency, 
suggesting that this increase in peak frequency is associated with 
dust in the photo-dissociation regions around the nebulae. The 
emissivity of spinning dust in these diffuse regions is of the same 
order as previous detections in the literature. Over the entire sky, 
the Commander solution finds more anomalous microwave emis¬ 
sion (AME) than the WMAP component maps, at the expense 
of synchrotron and free-free emission. Although the Commander 
model fits the data exceptionally well, as noted in Sect 9.4, the 
discrepancy is largely driven by differences in the assumed syn¬ 
chrotron spectrum and the more elaborate model of spinning 
dust designed to allow for the variation in peak frequency noted 
above. Future surveys, particularly at 5-20 GHz, will greatly im¬ 
prove the separation, since the predicted brightness between the 
two models disagrees substantially in that range. 

In polarization, synchrotron emission completely dominates 
on angular scales larger than 1° and frequencies up to 44 GHz. 
We combine Planck and WMAP data to make the highest signal- 
to-noise ratio map yet of the intensity of the all-sky polar¬ 
ized synchrotron emission at frequencies above a few giga¬ 
hertz, where Faraday rotation and depolarization are negligi¬ 
ble (Figs. 21 and Fig. 22). Most of the high-latitude polar¬ 
ized emission is associated with distinct large-scale loops and 
spurs, and we re-discuss their structure following the earlier 
study of Vidal et al. (2014) based on WMAP observations. We 
argue that nearly all the emission at -90° < / < 40° is part 
of the Loop I structure, and show that the emission extends 
much further in to the southern Galactic hemisphere than pre¬ 
viously recognized, giving Loop I an ovoid rather than circu- 
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lar outline. However, it does not continue as far as the ''Fermi 
bubble/microwave haze”, which probably rules out an associ¬ 
ation between the two structures. The South Polar Spur (SPS, 
see Fig. 21) is bordered by a polarized dust filament and as¬ 
sociated low-velocity Hi, analogous to the cold features long 
known to border Loop I around the North Polar Spur. We find 
two structures that could correspond to distant analogues of the 
radio loops, as predicted by Mertsch & Sarkar (2013), includ¬ 
ing one surrounding the Cygnus X star forming region, both of 
which are again associated with dust polarization. 

We identify a number of other faint features in the polar¬ 
ized sky, including a dearth of polarized synchrotron emission 
directly correlated with a narrow, roughly 20 ° long filament seen 
in Ha at high Galactic latitude, and also visible in the Faraday 
rotation map of Oppermann et al. (2012). Finally, we look for 
evidence of polarized AME; however, many AME regions are 
significantly contaminated by polarized synchrotron emission, 
and we find a 2cr upper limit of 1.6 % in the Perseus region. 


11.2. Polarized thermal dust emission 

Planck has produced the first all-sky map of the polarized 
emission from dust at submm wavelengths (Figs. 17 and 23). 
Compared with earlier ground-based and balloon-borne obser¬ 
vations (e.g., Benoit et al. 2004; Ward-Thompson et al. 2009; 
Matthews et al. 2009; Koch et al. 2010; Matthews et al. 2014) 
this survey is an immense step forward in sensitivity, coverage, 
and statistics. It provides new insight into the structure of the 
Galactic magnetic field and the properties of dust, as well as the 
first statistical characterization of one of the main foregrounds 
to CMB polarization. The wealth of information encoded in the 
all-sky maps of polarized intensity, F, polarization fraction, p, 
and polarization angle, presented in Planck Collaboration X 
(2015) is illustrated in Fig. 24. Here we summarize the main 
results from the data analysis by the Planck Consortium. The re¬ 
lease of the data to the science community at large will trigger 
many more studies. 


11.2.1. The dust polarization sky 

Planck Collaboration Int. XIX (2015) presents an overview of 
the polarized sky as seen by Planck at 353 GHz, the most sensi¬ 
tive Planck channel for polarized thermal dust emission, focus¬ 
ing on the statistics of p and ij/. At all Ah below 10^^ cm“^, p 
displays a large scatter. The maximum p, observed in regions of 
moderate hydrogen column density (Ah < 2 x 10 ^^ cm“^), is high 
( 7 >max ~ 20 %). There is a general decrease in p with increasing 
column density above Ah ~ 1 x 10^^ cm“^ and in particular a 
sharp drop above Ah ~ 10 ^^ cm“^. 

The spatial structure of ij/ is characterized using the angle 
dispersion function <S, the local dispersion of if/ introduced by 
Hildebrand et al. (2009). The polarization fraction is found to be 
anti-correlated with S. The polarization angle is ordered over 
extended areas of several square degrees. The ordered areas 
are separated by long, narrow structures of high S that high¬ 
light interfaces where the sky polarization changes abruptly. 
These structures have no clear counterpart in the map of the 
total intensity, /. They bear a morphological resemblance to 
features detected in gradient maps of radio polarized emission 
(lacobelli et al. 2014). 


11.2.2. The Galactic magnetic field 

The Planck maps of p and ij/ contain information on the mag¬ 
netic field structure. The data have been compared to syn¬ 
thetic polarized emission maps computed from simulations of 
anisotropic magnetohydrodynamical turbulence, assuming sim¬ 
ply a uniform intrinsic polarization fraction of dust grains 
(Planck Collaboration Int. XX 2015). The turbulent structure of 
the magnetic field is able to reproduce the main statistical prop¬ 
erties of p and ij/ that are observed directly in a variety of nearby 
clouds (dense cores excluded). The large-scale field orientation 
with respect to the line of sight plays a major role in the quanti¬ 
tative analysis of these statistical properties. This study suggests 
that the large scatter of p at Ah smaller than about 10 ^^ cm“^ is 
due mainly to fiuctuations in the magnetic field orientation along 
the line of sight, rather than to changes in grain shape and/or the 
efficiency of grain alignment. 

The formation of density structures in the interstellar 
medium involves turbulence, gas cooling, magnetic fields, and 
gravity. Polarization of thermal dust emission is well suited 
to studying the role of the magnetic field, because it images 
structure through an emission process that traces the mass of 
interstellar matter (Planck Collaboration XI 2014). The Planck 
I map shows elongated structures (filaments or ridges) that 
have counterparts in either the Stokes Q or U map, or in 
both, depending on the mean orientation. The correlation be¬ 
tween Stokes maps characterizes the relative orientation be¬ 
tween the ridges and the magnetic field. In the diffuse in¬ 
terstellar medium, the ridges are preferentially aligned with 
the magnetic field measured on the structures. This statisti¬ 
cal trend becomes more striking for decreasing column density 
and, as expected from the potential effects of projection, for in¬ 
creasing polarization fraction (Planck Collaboration Int. XXXII 
2014). Towards nearby molecular clouds the relative orientation 
changes progressively from preferentially parallel in areas with 
the lowest Ah to preferentially perpendicular in the areas with 
the highest Ah (Planck Collaboration Int. XXXV 2015). This 
change in relative orientation might be a signature of the for¬ 
mation of gravitationally-bound structures in the presence of a 
dynamically-important magnetic field. 

The relation between the structure of mat¬ 
ter and the magnetic field is also investigated in 
Planck Collaboration Int. XXXIII (2014), modelling the 
variations of the Stokes parameters across three filaments for 
different hypotheses on p. For these representative structures in 
molecular clouds the magnetic fields in the filaments and their 
background have an ordered component with a mean orientation 
inferred from Planck polarization data. However, the mean 
magnetic field in the filaments does not have the same orienta¬ 
tion as in the background, with a different configuration in all 
three cases examined. Planck Collaboration Int. XXXIV (2015) 
analyzes the magnetic field in a massive star forming region, the 
Rosette Nebula and parent molecular cloud, combining Faraday 
rotation measures from the ionized gas with dust polarized 
emission from the swept-up shell. This same methodology and 
modelling framework could be used to study the field structure 
in a sample of massive star forming regions. 

11.2.3. Dust polarization properties 

Galactic interstellar dust consists of components with different 
sizes and compositions and consequently different polarization 
properties. The relatively large grains that are in thermal equi¬ 
librium and emit the radiation seen by Planck in the submm also 
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Fig. 21. Synchrotron polarization amplitude map, P = ^jQ^ + at 30 GHz, smoothed to an angular resolution of 60', produced 
by a weighted sum of Planck and WMAP data as described in (Planck Collaboration XXV 2015). The traditional locii of radio 
loops I-IV are marked in black, a selection of the spurs identified by Vidal et al. (2014) in blue, the outline of the Fermi bubbles 
in magenta, and features discussed for the first time in (Planck Collaboration XXV 2015) in red. Our measured outline for Loop I 
departs substantially from the traditional small circle. 



Fig. 22. All-sky view of the angle of polarization at 30 GHz, rotated by 90° to indicate the direction of the Galactic magnetic field 
projected on the plane of the sky. The colours represent intensity, dominated at this frequency by synchrotron emission. The “drap¬ 
ery” pattern was obtained by applying the line integral convolution (LIC; Cabral & Leedom 1993) using an IDL implementation 
provided by Diego Falceta-Goncalves (http: //each. uspnet. usp. br/fgoncalves/pros/lic. pro). Where the field varies sig¬ 
nificantly along the line of sight, the orientation pattern is irregular and difficult to interpret. 
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Fig. 23. Dust polarization amplitude map, P = at 353 GHz, smoothed to an angular resolution of 10', produced by the 

diffuse component separation process described in (Planck Collaboration X 2015) using Planck and WMAP data. 



Fig. 24. All-sky view of the angle of polarization at 353 GHz, rotated by 90° to indicate the direction of the Galactic magnetic 
field projected on the plane of the sky. The colours represent intensity, dominated at this frequency by thermal dust emission. 
The “drapery” pattern was obtained by applying the line integral convolution (LIC; Cabral & Leedom 1993) using an IDL imple¬ 
mentation provided by Diego Falceta-Goncalves (http: //each. uspnet. usp. br/fgoncalves/pros/lic. pro). Where the field 
varies significantly along the line of sight, the orientation pattern is irregular and difficult to interpret. 
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extinguish and polarize starlight in the visible (Martin 2007). 
Comparison of polarized emission and starlight polarization on 
lines of sight probed by stars is therefore a unique opportunity to 
characterize the properties of polarizing grains. For this compar¬ 
ison, Planck Collaboration Int. XXI (2015) use P and 1 in the 
Planck 353 GHz channel and stellar polarization observations 
in the V band, the degree of polarization, pv, and the optical 
depth to the star, ry. Lines of sight through the diffuse interstellar 
medium are selected with comparable values of the column den¬ 
sity as estimated in the submm and visible and with polarization 
directions in emission and extinction that are close to orthogonal. 
Through correlations involving many lines of sight two ratios are 
determined, Rs/v = {Pll)l{pvl^v) and Rpjp = PIpv, the latter 
focusing directly on the polarization properties of the grains con¬ 
tributing to polarization. The first ratio, Rs/v^ is compatible with 
predictions based on a range of dust models that have been de¬ 
veloped for the diffuse interstellar medium (e.g., Martin 2007; 
Draine & Fraisse 2009). This estimate provides new empirical 
validation of many of the common underlying assumptions of 
the models, but is not very discriminating among them. The 
second ratio, Rp/p, is higher than model predictions by a factor 
of about 2.5. A comparable difference between data and model 
is observed for I/ry (Planck Collaboration Int. XXIX 2014). To 
address this, changes will be needed in the optical properties of 
the large dust grains contributing to the submm emission and 
polarization. 

The spectral dependence in the submm is also important 
for constraining dust models. In Planck Collaboration Int. XXII 
(2015) the Planck and WMAP data are combined to characterize 
the frequency dependence of emission that is spatially correlated 
with dust emission at 353 GHz, for both intensity and polariza¬ 
tion, in a consistent manner. At v > 100 GHz, the mean spec¬ 
tral energy distribution (SED) of the correlated emission is well 
fit by a modified blackbody spectrum for which the mean dust 
temperature of 19.6 K (derived from an SED fit of the dust total 
intensity up to 3000 GHz, i.e., 100yum) is adopted. It is found 
that the opacity has a spectral index 1.59 + 0.02 for polarization 
and 1.51 + 0.01 for intensity. The difference between the two 
spectral indices is small but significant. It might result from dif¬ 
ferences in polarization efficiency among different components 
of interstellar dust. Planck Collaboration Int. XXII (2015) also 
finds that the spectral energy distribution increases with decreas¬ 
ing frequency at v < 60 GHz, for both intensity and polarization. 
The rise of the polarization SED towards low frequency might be 
accounted for by a synchrotron component correlated with dust, 
with no need for any polarization of the anomalous microwave 
emission. 

11.2.4. Polarized dust and the CMB 

The polarized thermal emission from diffuse Galactic dust is 
the main foreground present in measurements of the polariza¬ 
tion of the CMB at frequencies above 100 GHz. The Planck sky 
coverage, spectral coverage from 100 to 353 GHz for HFI, and 
sensitivity are all important for component separation of the po¬ 
larization data. Planck Collaboration Int. XXX (2014) measures 
the polarized dust angular power spectra and over the 
multipole range 40 < ^ < 600 well away from the Galactic 
plane, providing cosmologists with a precise characterization of 
the dust foreground to CMB polarization. 

The polarization power spectra of the dust are well de¬ 
scribed by power laws in multipole, oc with exponents 
a = -2.42 + 0.02 for both the EE and BB spectra. The am¬ 
plitudes of the polarization power spectra are observed to scale 


with the average dust brightness as similar to the scal¬ 

ing found earlier for power spectra of 1 (Miville-Deschenes et al. 
2007). The frequency dependence of the power spectra for polar¬ 
ized thermal dust emission is consistent with that found for the 
modified blackbody emission in Planck Collaboration Int. XXII 
(2015). A systematic difference is discovered between the ampli¬ 
tudes of the Galactic B- and £’-modes, such that = 0.5. 

There is additional information coming from the dust TE and TB 
spectra. These general properties apply at intermediate and high 
Galactic latitude in regions with low dust column density. The 
data show that there are no windows in the sky where primordial 
CMB 5-mode polarization can be measured without subtraction 
of polarized dust emission. 


12. Summary and Conclusions 

This paper is an overview of the Planck 2015 release, summa¬ 
rizing the main features of the products being released and the 
main scientific conclusions that we draw from them at this time. 
Some of the highlights of this release are listed below. 

- Data from the entire mission are now used, including both 
temperature and polarization, and significant improvements 
have been made in the understanding of beams, pointing, cal¬ 
ibration, and systematic errors. As a result, the new products 
are less noisy, but even more importantly they are much bet¬ 
ter understood and the overall level of confidence is signifi¬ 
cantly increased. 

- The residual systematics in the Planck 2015 polarization 
maps have been dramatically reduced compared to 2013, 
by as much as two orders of magnitude in some cases. 
Nevertheless, on angular scales greater than 10°, systematic 
errors in the polarization maps between 100 and 217 GHz 
are still non-negligible compared to the expected cosmolog¬ 
ical signal. It was not possible, for this data release, to fully 
characterize the large-scale residuals due to these system¬ 
atic errors from the data or from simulations. Therefore all 
results published by the Planck Collaboration in 2015 have 
used CMB polarization maps that have been high-pass fil¬ 
tered to remove the large angular scales. Users of the Planck 
CMB maps are warned that they are not useable for cosmo¬ 
logical analysis at Z’i30. 

- A large set of simulations accompanies the release, includ¬ 
ing up to 10000 realizations of signal and noise; this has 
been used to test and verify methods of analysis, and also to 
estimate uncertainties. 

- We measure the amplitude and direction of the Solar dipole 
to the best precision so far. 

- One of the most notable improvements in this data set is that 
now LFI, HFI, and WMAP agree to within a few tenths of 
a percent on angular scales from the dipole through the first 
acoustic peak. 

- Polarization is a new product in this release, and especially 
on large angular scales systematic effects are not yet fully 
controlled. This is an area where we expect to make signifi¬ 
cant improvements in the coming months. Nevertheless, our 
polarization data are already making important contributions 
in a variety of analyses. 

- More specifically, we are able to use the TE and EE angular 
power spectra at small scales (and to a more limited extent, at 
large angular scales) over the full sky, reaching the expected 
sensitivity. This allows us to estimate cosmological parame¬ 
ters independently of TT, and in combination with TT. 
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- At large angular scales, we are now able to use Planck-ovi[y 
products to carry out cosmological analysis. Specifically, we 
can estimate the optical depth of reionization, r, indepen¬ 
dently of other experiments. The value of r is smaller than 
found in previous determination, implying later reionization. 

- Foregrounds can be separated effectively over larger areas of 
the sky, allowing more sky to be used for cosmology, and 
producing high-quality maps of synchrotron, free-free, spin¬ 
ning dust, thermal dust, and CO emission. 

- Our 2015 results for cosmology are very consistent with our 
2013 results, but with smaller uncertainties, and covering a 
greater range of science. 

- Our best-fit 2015 cosmological parameters comfirm the basic 
6-parameter ACDM scenario that we determined in 2013. 
There is no compelling evidence for any extensions to the 
6-parameter model, or any need for new physics. Depending 
somewhat on what data combinations are used, five of the six 
parameters are now measured to better than 1 % precision. 
Areas that were in “tension” in 2013 (erg and weak galaxy 
lensing), have been confirmed to remain in tension today, 
although the disagreement is lessened when only particular 
subsets of the external data are considered. 

- Using only Planck data, we find that the Universe is fiat to 
0.7 % (1 cr). Including BAO data, the constraint tightens to a 
remarkable 0.25 %. 

- Using the Planck temperature data over the whole sky, 
together with our recent work combining Planck and 
BICEP2/Keck data, we have obtained the best current upper 
limits on the tensor-to-scalar ratio obtained to date. 

- We have obtained improved limits on primordial non- 
gaussianity (/nl), which are about 30% tighter before, 
reaching the expected sensitivity of Planck when including 
polarization. 

- Models of inflation are more tightly constrained than ever 
before, with the simplest 0" models being ruled out for n>2. 

- We have obtained the tightest limits yet on the amplitude of 
primordial magnetic fields. 

- Planck's measurement of lensing of the CMB has the highest 
signal-to-noise ratio yet achieved, 40 cr. 

- The second Planck catalogues of compact sources, Sunyaev- 
Zeldovich clusters, and Galactic cold clumps, are larger than 
the previous ones and better-characterized in terms of com¬ 
pleteness and reliability. 

Planck continues to provide a rich harvest of data for cos¬ 
mology and astrophysics. 
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